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ABSTRACT 


Deactivation  coefficients,  yw  of  vibrationally-excited  N^X^.v)  on  alloys  of  aluminum, 
stainless  steel  alloys,  a  titanium  alloy,  gold,  Pyrex  glass  and  Teflon  were  determined  from 
measurements  of  the  wall  relaxation  rate  coefficients  for  the  v  =  1  to  4  or  5  hot  bands  using 
CARS  A  water-cooled  positive  column  discharge  was  used  to  populate  the  vibrational 
states.  The  nitrogen  gas  flow  rate  was  adjusted  to  optimize  the  observed  decay  of  the 
populations  The  excited  gas  was  exposed  to  the  tubular  specimen  surface  for  precisely 
controlled  residence  times.  The  BOXCARS  configuration  was  used  to  measure  only  the  gas 
that  was  exposed  to  the  specimen  surface.  Extensive  measurements  on  Pyrex  gave  y,  values 
between  2  .4*  lO4  and  6.7*  1CT4  depending  on  the  treatment  history  of  the  surface.  The  values 
for  y4  ranged  from  2.9*  10"4  for  the  AMS  4943D  alloy  of  titanium  to  approximately  unity  for 
the  AMS  3 12  stainless  steel  alloy.  The  low  value  for  titanium  can  be  attributed  to  the  oxide 
layer.  The  variation  of  k^,  with  v  was  linear  or  nearly  linear  in  all  cases  with  slopes  lower  in 
most  cases  than  the  rate  of  increase  of  the  vibrational-translational  exchange  rate  with  v. 
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INTRODUCTION 

Vibrationally  excited  molecules  play  various  roles  in  gas  kinetics  and  plasma 
chemistry.  Of  the  seven  diatomic,  homonuclear  molecules  occurring  at  normal  temperatures, 
nitrogen  holds  interest  due  to  its  importance  in  combustion,  biology,  microelectronics,  high 
speed  and  high  altitude  aircraft  and  spacecraft  reentry.  Vibrational  energy  in  the  ground 
electronic  state  of  nitrogen  is  important  in  the  short-lived  pink  afterglow  production  where 
Penning  ionization  and  the  rapid  e-V  exchange  leads  to  considerable  heat  loss  and  N2(B2IU+) 
pumping.1  Vibrational  energy  is  important  in  promoting  plasma  chemistry  where  it  is  an 
available  energy  source  for  molecular  dissociation.  Here  radicals  or  atoms  are  delivered 
where  needed  as  in  plasma-enhanced  chemical  vapor  deposition  (PECVD) 2  Published 
research  into  the  catalytic  chemistry  of  active  nitrogen  in  low  pressure  post-discharge  gas 
revealed  the  presence  of  a  large  concentration  of  excited  ground  state  nitrogen  N2(X]Ig+,  v) 
accompanying  N(4S3/2)  atoms 3  Accurate  calorimetric  measurements  of  atom  recombination 
require  correction  for  the  vibrational  energy  deposition.  Accurate  theoretical  models  of 
vibrational  distribution  require  specific  state-selective  measurements,  which  have  been 
developed4-5 

Vibrational  energy  transfer  ofN^X’S^,  v)  to  solid  surfaces  is  addressed  in  this  work. 
Population  distributions  and  loss  rates  in  the  presence  of  a  surface  were  studied.  These 
results  extend  the  knowledge  about  diatomic  molecular  wall  loss  rates  to  higher  v-states  than 
is  currently  available.  This  information  can  influence  the  design  of  gas  lasers  and  improve  our 
understanding  of  all  reactive  gaseous  environments,  such  as  in  combustors,  in  plasma 
processing  chambers  for  microelectronics  and  on  high  temperature  exterior  surfaces  of  high 
speed  aircraft  In  this  work,  vibrationally  excited  nitrogen  was  exposed  to  tubular  reactor 
surfaces  of  common  construction  materials.  This  list  includes  neat  Pyrex,  Pyrex  exposed  to 
active  nitrogen,  type  3003  aluminum,  types  304  and  AMS  312  stainless  steel,  AMS  4943D 
titanium,  99.999%  pure  gold  and  Teflon.  The  tubular  reactor  elements  had  an  internal 
diameter  of  about  1.5  cm  and  were  cut  into  specimens  about  12  cm  long. 
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A  flowing  afterglow  created  in  a  water-cooled,  dc  electric  discharge  reacted  with  each 
surface  for  precisely-stepped  residence  times.  The  glow  discharge  provided  a  low 
translational  temperature  and  a  high  vibrational  excitation  with  insignificant  atom  and  charged 
species  densities.  This  gave  spectra  which  had  isolated  vibrational  bands  with  significant 
rotational  line  intensities  limited  to  approximately  J  =  21.  The  vibrational  wall  relaxation  rates 
were  measured  for  the  v  =  1  to  either  the  v  =  4  or  5  vibrational  “hot”  bands  by  probing  the 
axial  flow  with  a  small-volume  optical  probe  using  coherent  anti-Stokes  Raman  spectroscopy 
(CARS).  Although  signals  from  the  v  =  6  band  were  observed,  useful  measurements  above 
the  v  =  5  band  were  virtually  impossible  under  our  conditions.  In  essence  this  means  those 
populations  were  too  small  to  have  an  impact  on  the  experiment  as  a  whole.  Nonequilibrium 
distributions  were  characterized  by  rotational  temperatures  of  -350  K  and  vibrational 
temperatures  of  -2000  K.  Numerous  experiments  were  made  on  some  surfaces  to  determine 
repeatability.  At' the  same  time,  a  slow  change  in  the  surface  due  to  the  afterglow  was 
observed  for  every  surface.  Measurements  were  made  on  the  axis  of  the  gas  flow  and 
referred  to  the  surface  by  calculation. 

Observed  relaxation  rate  k^,  which  include  both  homogeneous  and  heterogeneous 
deactivation,  and  the  “apparent”  deactivation  coefficients  yva  for  the  surfaces  are  reported  for 
each  measured  v  -►  v  +  1  band.  The  vibrational  distribution  in  the  bulk  of  the  gas  depends 
on  the  relaxation  rates  and  radial  diffusion.  Relaxed  molecules  must  diffuse  to  the  probe 
volume  before  the  effect  can  be  seen.  Moreover,  the  observed  deactivation  due  to  the  surface 
is  significantly  modified  by  the  fast  vibrational-vibrational  (V-V)  exchange  processes  via 
diffusion.  These  processes  soften  the  effect  of  wall  relaxation  by  the  redistribution  of  the 
populations  up  the  vibrational  energy  level  ladder.  This  process  is  strong  when  there  is 
significant  populations  in  the  low-lying  vibrational  levels  such  as  produced  by  a  gaseous 
discharge.  For  the  low  lying  v-levels  studied  here  (<5),  the  vibrational-translational  (V-T) 
exchange  rates  are  insignificant.  In  order  to  extract  the  “true”  wall  relaxation  rates  k^, 
from  ky,  a  Master  Equation  Model  of  N2  was  fitted  to  the  data  for  four  specimens  with  the 
kv<waii)  as  fitting  parameters.  This  model  includes  V-V  and  V-T  energy  exchange  processes 
as  well  as  wall  exchange  interactions.  It  uses  the  Schwartz-Slawsky-Herzfeld  (SSH)  law  of 
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scaling  relaxation  rates  and  the  development  of  the  model  closely  follows  that  of  Bailey6  and 
Capitelli 7  From  this  procedure,  “true”  values  of  and  yv  were  determined.  As  described 
in  a  later  section,  this  fitting  effort  was  greatly  complicated  by  the  fact  that  the  model 
computed  populations  out  to  v  =  40  while  data  was  available  only  up  to  v  =  4  or  5.  This 
required  an  extrapolation  procedure  to  constrain  the  fitting  of  the  populations  above  v  =  4  or 
5.  It  was  found  that  many  of  the  data  sets  could  not  be  adequately  constrained  to  permit 
unambiguous  fittings.  Nevertheless,  it  is  suggested  that  the  four  successful  fittings  provide 

information  that  can  be  used  to  provide  reasonable  corrections  of  the  yVB  for  all  the  specimens 
to  viable  estimates  of  the  “true”  yv. 
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EXPERIMENTAL  DETAILS 


In  an  experimental  apparatus  similar  to  that  used  by  Black  et  al.,8  the  vibrationally- 
excited  molecules,  generated  in  the  low  pressure  dc  discharge,  were  carried  in  the  afterglow 
through  a  15-cm-long  tubular  reactor  in  which  diffusion  carried  them  to  the  surface.  Black 
et  al.  used  mostly  a  thermal  source  and  a  100-cm-long  reactor.  Because  of  this  long  flow  tube 
length,  they  had  to  correct  for  Poiseuille  flow.  Moreover,  they  could  only  achieve  sufficient 
excitation  to  observe  the  v  =  1  -*  0  transition  using  spontaneous  anti-Stokes  Raman 

scatttering.  A  fraction  yv  of  the  excited  molecules  which  collide  with  a  surface  having  a 
catalytic  capability  result  in  the  transfer  of  vibrational  energy  to  the  surface  vibrational 
modes.9  The  population  distribution  remaining  in  the  gas  was  measured  after  stepped 
residence  times  of  contact  of  the  gas  with  the  surface.  The  pressure  chosen  for  these 
measurements  was  typically  between  10  to  20  Torr.  This  pressure  was  sufficient  to  provide 
adequate  CARS  signals  and  low  enough  to  permit  efficient  steady  state  production  of 
vibrationally-excited  nitrogen  in  the  discharge  From  the  observed  decay  of  the  vibrational 
populations  with  time,  the  total  rate  coefficients  ky  were  determined  from  which  the  apparent 

wall  deactivation  efficiencies  yva  were  calculated  for  each  vibrational  level.  The  total  rate 
coefficients  are  independent  of  the  size  of  the  probe  volume  as  well  as  the  radial  location  in 
the  tubular  specimen.  They  do  depend,  however,  on  the  combined  effects  of  collisions  with 
the  surface  and  the  gas-gas  collisions  For  the  conditions  of  these  studies,  the  latter  is 

dominated  by  V-V  processes.  For  given  values  of  “true”  yv,  the  observed  rate  coefficients 

ky  depend  on  the  temperature,  pressure  and  tube  radius.  The  yv,  however,  are  characteristics 
of  the  surface  which  are  determined  by  surface  activity,  temperature,  the  colliding  molecule, 
and  possibly  pressure.  By  fitting  the  Master  Equation  Model  to  these  population  data, 
extraction  of  k^^,  free  of  homogeneous  effects  was  attempted.  However,  because  the 
coupled  equations  are  stiff  and  the  fitting  process  is  highly  dependent  on  the  fact  that  the  rate 
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of  change  of  population  in  one  level  is  highly  coupled  to  neighboring  populations  via  the  V- V 
processes,  this  effort  gave  mixed  results. 

A  flowing  gas  handling  system  carried  nitrogen  (Air  Products  ultra-pure  carrier  grade 
99.999%  pure)  through  a  water-cooled  dc  electric  discharge  operating  at  about  2  Td  with  40 
cm  of  the  positive  column  exposed  to  the  flow.  The  flowing  afterglow  passed  around  a  right- 
angle  bend  and  then  through  the  tubular  reactor  specimen.  Tylan  General  FC-260  flow 
controllers  maintained  steady  mass  flow  to  in  the  range  of  0  to  2000  seem  calibrated  to  within 
0.2%.  The  gas  supply  system  was  all  stainless  steel  except  for  the  Pyrex  discharge  tube,  the 
source  tube  (which  carried  the  active  nitrogen  to  the  specimen  tube)  and  a  0.8  m  section  of 
flexible  polyethylene  tubing  at  the  input  of  the  discharge  tube  which  prevented  the  plasma 
from  striking  up-stream.  All  experiments  were  accomplished  between  0.5  and  23  Torr  (67 
to  3066  Pa)  with  most  measurements  at  14  Torr.  In  Fig.l,  the  positive  column  discharge 
pumps  the  vibrational  states  to  nonequilibrium.  The  vibrational  temperature  of  the  gas 
entering  the  specimen  tube,  T„  ranged  from  about  1800  to  2800  K  while  the  gas  temperature, 
Tg,  was  between  about  3 15  to  380  K.  Either  a  dc  discharge  or  a  cw  microwave  source  could 
have  been  used  to  preferentially  pump  N2(X,v).10  The  dc  discharge  avoided  high 
concentrations  of  N(4S)  atoms  and  N2  B-  and  A-state  molecules.  Recombination  of  the 
residual  atom  concentration  produced  the  weak  straw-colored  Lewis-Rayleigh  afterglow.111 
A  spectroscopic  study  of  this  glow  with  residence  time  provided  a  measure  of  atom 
recombination  on  the  wall  and  in  the  gas  and  gave  a  rough  estimate  of  N  density.12  The  mole 
fraction  of  N,  N2(B  and  A)  species  were  subsequently  found  to  be  about  0.0005.  This 
was  less  than  the  theoretical  density  of  v  =  9  or  10.  Thus,  molecules  not  in  the  electronic 
ground  state  and  atoms  had  only  a  negligible  effect  on  the  results  reported  in  these  studies. 
The  transit  time  of  the  gas  in  the  source  tube  was  no  less  than  20  ms.  Because  of  this,  all 
charged  species  had  sufficient  time  to  become  neutralized  well  before  they  came  in  contact 
with  the  specimen  surface. 

Exiting  the  discharge,  the  gas  turned  90  degrees,  past  a  Wood's  horn,  into  the  source 
tube  where  a  pink  afterglow  terminated  The  line-of-sight  from  the  discharge  tube  terminated 
into  the  Wood’s  horn  in  order  to  prevent  near-UV  emissions  from  the  plasma  exciting  the  gas 
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Figure  1 .  Schematic  diagram  of  discharge  tube  and  the  main  probe  chamber 
for  measuring  the  relaxation  rates  as  a  function  of  Z. 
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entering  the  reactor  specimen  tube.  During  the  transit  to  the  specimen  tube,  the  gas  cooled 
and  attained  steady  state  conditions  giving  a  steady  population  at  the  exit  of  the  source  tube 
The  position  of  the  source  tube,  which  was  sealed  through  an  O-ring  gland  in  the  top  of  the 
optically-accessible  main  chamber,  could  be  positioned  to  ±0. 1  mm,  providing  very  accurate 
control  of  the  residence  time  of  the  gas  in  the  specimen  tube.  Brewster-angle  windows  were 
mounted  on  15-cm-long  tubular  extensions  on  the  chamber  which  permitted  focusing  of  the 
high  peak-power  CARS  beams  in  the  chamber  without  damaging  the  windows.  There  was 
also  a  pair  of  side  windows  for  emission  experiments.  The  CARS  signals  were  generated  by 
probing  the  gas  on  the  axis  of  the  specimen  tube  in  a  1-cm  gap  between  bottom  of  the 
specimen  and  the  top  of  a  catch  tube.  The  catch  tube  was  connected  to  a  vacuum  pump 
through  a  set  of  three  valves  arranged  to  provide  for  on/off,  coarse  and  fine  control  of  the 
pumping  rate  on  the  main  chamber.  Pressure  was  regulated  by  manual  butterfly  and  needle 
valves  some  1 50  cm  distance  after  the  specimen  tube  and  monitored  with  a  pair  of  Wallace 
and  Tieman  gauges  (0  to  20  and  0  to  200  Torr  absolute).  Additional  6-mm-diameter  vents 
in  the  floor  of  the  main  chamber  constrained  recirculation  of  the  gas.  A  computer-generated 
analog  signal  controlled  the  mass  flow  controller.  Thermocouple  and  Convectron  gauges 
were  incorporated  for  base  pressure  and  helium  leak  checking  of  the  system,  respectively. 

At  more  than  a  meter  upstream  of  the  discharge  tube,  a  mixing  chamber  that  had  three 
inlet  ports  was  used  to  mix  nitrogen  with  argon  or  helium  as  needed.  The  system  was 
returned  to  atmospheric  pressure  with  pure  nitrogen  or  argon.  Before  each  experiment,  each 
specimen  was  degreased  with  methylene  chloride  in  a  sonic  cleaner  before  insertion  in  the 
main  chamber.  The  source  and  specimen  surfaces  were  cleaned  with  an  argon  discharge 
around  2  Torr  at  60  mA.  The  pressure  was  varied  between  0.25  and  10  Torr  during  this 
procedure  to  promote  out-gassing  of  any  trapped  gases.  No  oxygen  or  water  could  be 
detected  in  gas  samples  taken  from  the  system  and  analyzed  locally  by  Fourier  transform  mass 
spectrometry.13  The  supply  and  exhaust  systems  were  periodically  baked  out  assisted  by  a 
combination  of  heat  tapes  and  heat  gun  with  either  flowing  gas  or  pumping  at  the  1  pm  base 
pressure  for  six  to  eight  hours.  Base  pressure  was  accomplished  with  an  isolated  secondary 
system  including  an  Edwards  mechanical  pump  fitted  with  a  molecular  sieve  at  the  inlet  port. 
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The  primary  flow  system  was  pumped  by  a  Welch  mechanical  pump  fitted  with  an  oil  mist 
filter  and  vented  outside  the  lab.  At  no  time  were  primary  and  secondary  pumps  operated 
together.  Moreover,  the  primary  pump  was  never  opened  at  or  below  its  base  pressure  of  25 
pm  and  was  always  operated  with  flowing  gas. 

Mass  flow  was  regulated  to  a  fixed  value  for  each  experiment  between  600  and  1200 
seem  ±0.2%.  This  provided  an  axial  transport  speed  between  about  6  and  1 1  m/s  in  the 
specimen  tube.  This  gave  from  12  to  30  ms  maximum  residence  time  of  the  vibrationally 
active  gas  in  the  specimen  tube  and  more  than  20  ms  in  the  source  tube.  Because  the 
specimen  tube  was  short,  plug  flow  was  assumed,  thereby  requiring  no  correction  for 
Poiseuille  flow. 

CARS  is  an  effective  probe  of  ground  electronic  state  vibrational  levels  in  homopolar 
diatomic  molecules.  Lack  of  an  absorption  spectrum  due  to  the  absence  of  electric  dipole 
transitions  compels  the  use  of  a  Raman  probe.  At  low  pressures  the  nonlinear  process 
produces  little  signal  generation  advantage  over  linear  Raman  spectroscopy;  however,  the 
nonlinear  process  generates  a  signal  that  is  confined  to  a  coherent  beam  which  greatly 
simplifies  the  signal  collection.  On  the  other  hand,  two  complications  demand  careful 
attention  The  change  in  the  third  order,  nonlinear  susceptibility  in  the  CARS  process  is 
proportional  to  the  difference  between  the  populations  of  the  two  quantum  levels  probed  by 
the  laser  beams.  This  requires  an  analysis  step,  discussed  in  the  Analysis  section,  that  extracts 
fractional  populations  from  signal  data  containing  population  differences.  Also,  stimulated 
Raman  scattering  (SRS)  simultaneously  pumps  population  from  the  lower  to  the  upper  level, 
which  introduces  a  laser-power-dependent  systematic  shift  in  the  signal  measurements.  This 
“saturation”  effect  has  been  the  subject  of  an  investigation  similar  to  studies  reported  here.14 
However,  it  was  necessary  to  independently  determine  the  correction  for  the  intermediate 
peak  laser  power  levels  used  in  these  studies.15  This  correction  requires  a  calibration  step 
which  is  briefly  described  in  the  Analysis  section. 

The  pulsed  lasers  used  in  the  CARS  system  consisted  of  an  injection  seeded,  doubled 
Spectra-Physics  DCR-3  Nd  YAG  laser  which  pumped  a  Lumonics  HyperDye  300  tunable  dye 
laser.  The  YAG  laser  produced  a  7  ns  pulse  and  provided  a  CARS  pump  beam  with  a  single 
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longitudinal  mode  linewidth  of  0.006  cm'1  (0.00018  nm,  180  MHZ).  The  dye  laser,  using 
sulforhodamine  640dye,  produced  the  5  ns,  multimode  CARS  probe  beam  with  a  linewidth 
ofO.l  cm'^O.OCB?  nm,  2  GHz)  over  the  range  of  16638  to  16393  cm'1  (601  to  610  nm).  It 
was  found  that  22.33  mg  dye  per  liter  methanol  dilution  provided  a  slowly  varying  power 
curve  over  this  range  Figure  2  is  a  diagram  of  the  optical  system  used  in  this  work.  The 
figure  shows  how  the  appropriate  timing  delay  of  the  YAG  beams  was  accomplished  through 
a  series  of  prisms.  The  YAG  and  dye  laser  beam  energies  were  adjusted  using  half-wave 
plates  in  combination  with  polarizing  prisms.16  The  main  pump  beam  was  split  into  two  equal 
energy,  parallel  beams  about  2  cm  apart  with  a  beam-splitter  mirror  and  aligned  symmetrically 
on  the  focusing  lens.  The  probe  beam  was  increased  in  diameter  through  a  beam  expanding 
telescope  in  order  to  approximately  match  the  dye  beam  waist  to  the  waist  sizes  of  the  two 
pump  beams  in  the  focal  volume  of  the  lens.  The  probe  beam  was  aligned  in  a  plane 
perpendicular  to  the  plane  of  the  pump  beams  so  as  to  satisfy  the  BOXCARS  phase  matching 
condition.17  The  pump  beams  were  elliptical  in  cross  section  with  major  and  minor  axis 
dimensions  at  full  width  at  half  maximum  (FWHM)  of 300  pm  and  100  |im  respectively.  The 
waist  diameter  of  the  dye  beam  was  70  pm  FWHM.  These  dimensions  were  determined  by 
knife-edge  test  and  bum  patterns.  The  length  of  the  overlap  or  probe  volume  of  the  three 
beams  was  measured  to  be  about  5  mm  by  translating  a  thin  argon  sheet  through  the  volume 
while  measuring  the  air  signal  The  tolerance  on  the  positioning  of  the  probe  volume  was  ±2 
to  3  mm  due  to  the  length  of  the  probe  volume  and  the  15-mm  tube  diameter.  This 
positioning  was  accomplished  by  translating  the  main  chamber  parallel  to  the  probe  beams 
while  measuring  the  CARS  signal  for  the  v  =  1  to  2  transition.  Beam  alignment  was  checked 
before  each  experiment  by  peaking  the  CARS  signal  with  the  dye  beam  focusing  lens.  During 
chamber  maintenance,  beam  crossing  alignment  was  achieved  by  passing  the  three  attenuated 
beams  through  a  pin  hole  The  BOXCARS  beam  arrangement  was  chosen  because  it 
permitted  precise  location  of  the  measurement  volume.  A  scheme  using  60°  rotated 
polarization  in  each  beam,  available  for  suppressing  the  nonresonant  background,  was  not 
used  because  the  signal  is  reduced  by  one  third.18  Since  the  nonresonant  background  in  these 
measurements  was  very  weak,  this  signal  reduction  was  deemed  too  great  considering  the 
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diminutive  size  of  the  signals  from  the  low  pressure  gas.  All  beams  were  horizontally 
polarized  parallel  with  the  optical  table. 

Beam  powers  incident  on  the  main  chamber  were  tracked  with  PIN  diodes  calibrated 
between  each  run  within  each  experiment.  Diode  counts  linearly  tracked  the  incident  fluence 
and  systematic  energy  deviations  in  the  laser  pulses.  While  the  diode  signal  counts  were 
proportional  to  the  incident  pulse  energy,  average  power  measurements  were  more 
convenient  to  use  because  the  average  diode  signal  was  used  in  normalizing  the  average 
CARS  signal.  These  diodes  were  an  important  nexus  between  saturation  calibration  and  the 
deactivation  experiments.  The  count-to-mW  ratios  of  each  deactivation  experiment  and  the 
associated  saturation  calibration  experiment  were  used  to  convert  the  power  in  the 
deactivation  measurements  to  the  equivalent  power  used  in  the  saturation  measurements.  The 
degree  of  saturation  S  was  then  calculated  and  returned  to  correct  the  raw  CARS  signal 
values. 

Instead  of  using  a  high  pressure  argon  gas  cell  to  reference  the  measurements,  which 
was  really  impractical  for  our  system  besides  being  notoriously  difficult  to  keep  aligned,  the 
averages  of  the  outputs  of  the  laser  PIN  diodes  were  used  to  normalize  the  CARS  signals. 
These  integrated  signals  tracked  power  changes  in  the  two  beams  within  a  run  as  well  as 
between  runs  within  an  experiment.  During  any  given  deactivation  or  saturation  run,  the 
standard  deviation  of  the  mean  of  the  output  of  the  YAG  laser  over  200  pulses  was  well 
under  1  %.  Due  to  mode  beating  in  the  dye  laser,  there  was  no  significant  correlation  between 
the  fluctuations  in  the  dye  laser  output  and  the  YAG  laser  output  Thus,  assuming  that  there 
was  no  significant  correlation  between  the  susceptibility  fluctuations  and  the  dye  laser 
fluctuations,  the  normalization  of  the  CARS  signal  was  accomplished  by  dividing  the  average 
CARS  signal  by  the  individual  averages  of  the  laser  diode  signals.  The  normalization  and 
saturation  correction  are  detailed  in  the  Analysis  section.  In  addition,  signal  changes  due  to 
the  dye  spectral  power  curve  were  removed  initially  by  calibrating  the  wavelength  response 
of  the  dye  laser  PIN  diode.  It  was  found  subsequently  that  the  combined  spectral  response 
of  the  receiver  optics,  spectrometer  and  photomultiplier  tube  (PMT)  was  flat  over  the 
operating  range  and  needed  no  correction. 


11 


The  receiver  optics  consisted  primarily  of  a  500-mm-focal-length  planoconvex  lens, 
a  right-angle  prism  for  spatial  separation  of  signal  and  laser  beams,  two  p-polarization 
dichroic  mirrors,  two  p-polarization  dichroic  filters,  and  a  39-mm,  best-form  lens  for  imaging 
the  CARS  beam  on  the  entrance  slit  of  a  1-m  Jarrell-Ash  grating  spectrometer.  The 
spectrometer  was  used  primarily  as  an  additional  filter  for  rejection  of  any  pump  beam  leakage 
light.  An  EMI  9816KB  PMT  performed  CARS  signal  detection  typically  operating  at  2100 
V.  The  rejected  beams  were  directed  to  beam  dumps  except  for  about  4%  reflected  from  the 
hypotenuse  of  two  prisms  toward  the  PIN  diodes. 

Data  acquisition  and  dye  laser  tuning  was  computer  controlled  through  a  Keithley 
Data  Acquisition  System  Model  570.  The  data  acquisition  and  control  (DAqCt)  software  was 
written  in  GW  BASIC  which  incorporated  Keithley  Soft500  BASIC  libraries  for  controlling 
the  Keithley  570.  The  final  DAqCt  program  version,  named  SB9.BAS,  consisted  of  three 
separate  programs:  SB9,  SB9KL1  and  SB9KL2.  This  software  had  several  features 
including  the  capability  for  manually  controlling  the  dye  laser  and  the  flow  controllers, 
automatically  tuning  the  dye  laser  from  a  script  file  and  carrying  out  the  acquisition  of  a 
continuous  spectrum  or  a  skip  spectrum  from  a  script  file. 

Gated  detection  of  the  PMT  and  PIN  diode  signal  pulses  was  used  to  minimize 
background  signals  and  noise  from  ambient  light  and  from  radio  frequency  interference 
generated  primarily  in  the  YAG  flash  lamp  power  supply.  The  system  trigger  pulse  was 
generated  with  a  PIN  diode  located  behind  the  beam  splitter  for  the  dye  laser.  This  trigger 
pulse  was  sent  to  a  LeCroy  Model  821  Discriminator  to  generate  a  NIM  standard  pulse  which 
triggered  a  Hewlett-Packard  (H-P)  Model  3012B  pulse  generator.  The  output  from  this 
generator  was  attenuated  and  split  into  a  clock  start  pulse  and  three  -1  6-V  gate  pulses,  each 
70  ns  long.  The  gate  pulses  provided  the  gating  of  the  LeCroy  Model  2249SG  fast  pulse 
analog-to-digjtal  converter  (ADC).  The  ADC  digitized  the  integrated  signal  pulses  from  the 
PMT  (i.e.,  the  CARS  signal)  and  from  the  YAG  laser  and  dye  laser  PIN  diodes.  These  pulses 
were  appropriately  delayed  with  coax  delay  cables  such  that  they  arrived  at  the  ADC  within 
the  gate  pulses.  Each  pulse  was  passed  first  through  a  10-step  Phillips  Scientific  Model  771 
pulse  amplifier.  The  gain  was  set  to  maximize  the  utilization  of  the  ADC’s  1000-count 
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dynamic  range  consistent  with  the  fluctuations  of  the  pulse.  The  CARS  signal  was 
additionally  attenuated  using  a  motor-controlled  neutral  density  filter  wheel  mounted  in  front 
of  the  PMT.  This  was  necessary  because  of  the  more  than  1000  to  1  range  of  the  CARS 
signals.  Both  the  amplifiers  and  the  neutral  density  filters  were  calibrated  to  accuracies  better 
than  0. 1%.  The  YAG  laser  pulse  noise  level  of  1%  relative  standard  deviation  was  the  lowest 
of  the  three  signals.  The  dye  laser  presented  about  4%  to  6%  noise  while  the  noise  on  the 
CARS  signal  was  typically  between  22  to  33%  relative  standard  deviation.  Each  laser  PIN 
diode  had  an  adjustable  optical  attenuator  that  greatly  alleviated  spatial  alignment 
requirements  and  removed  channel  spectra  in  the  dye  diode.  The  attenuators  consisted  of  a 
pair  of  telescoping  brass  tubes  with  an  opal  glass  window  at  the  interior  end  of  each  tube. 

A  fifth  pulse  from  the  first  H-P  pulse  generator  triggered  a  second  H-P  3012B  pulse 
generator  which  produced  a  2.3  ms  pulse  for  gating  the  Keithley  570.  The  computer  reads 
the  ADC  after  each  laser  pulse  as  many  times  as  was  programmed  for  that  deactivation  data 
point.  After  each  data  point,  the  program  totaled  the  data  and  stored  the  value  in  a  file.  The 
wavelength  was  stepped  and  the  process  repeated.  Later,  the  average  signal  values  and 
uncertainties  were  then  calculated 

The  water  cooling  of  the  nitrogen  glow  discharge  kept  the  gas  temperature  at  the 
specimen  tube  below  400  K.  This  discharge  tube  was  connected  to  a  manually-operated,  high 
voltage  dc  power  supply  through  a  high-power  50  kQ  ballast  resistor.  Typical  voltages  were 
about  8.8  kV  at  the  power  supply  giving  about  4.3  kV  at  the  discharge  terminals.  The 
discharge  current  was  set  to  90  mA  and  monitored  periodically  throughout  the  experiment. 
Hollow  Kovar  electrodes  were  used  in  the  discharge  tube. 
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THEORY 


CARS  Probe  Technique 

The  CARS  technique  was  originally  developed  for  the  combustion  community. 
Combustion  is  often  accompanied  by  fluorescence  emissions,  high  luminous  emissions  due 
to  soot  formation,  high  temperatures,  chemical  reactivity,  turbulence  and  sometimes  high 
transport  velocities.  Physical  probes  can  disturb  the  chemistry  or  the  population  distributions, 
they  are  not  generally  state  selective  and  they  may  be  unable  to  physically  survive  the 
environment.  An  optical  probe  can  overcome  many  of  these  difficulties.  CARS  has  the 
highly-desired  property  of  producing  a  large,  coherent  optical  signal  beam  on  the  short 
wavelength  side  (i.e.,  blue  side)  of  the  incident  laser  wavelength  thereby  removing  the  signal 
from  the  strong,  potentially  interfering  fluorescence  and  incandescence  emissions  that  are  on 
the  red  side  of  the  laser  wavelength.  These  emissions  are  typically  generated  when  a  laser 
beam  is  incident  on  a  complex  combustion  medium.  This  is  in  contrast  to  the  spontaneous 
Raman  or  laser-induced  fluorescence  techniques  wherein  the  signals  are  usually  on  the  red 
side  of  the  laser  wavelength.  These  techniques  also  have  the  disadvantage  of  requiring  signal 
detection  over  a  significant  solid  angle  which  introduces  spatial  resolution  restrictions  and 
limitations  of  use  on  highly  luminous  media.  CARS  has  the  disadvantages  of  requiring  two 
optical  ports  that  can  withstand  high  peak  power,  short  optical  pulses  and  a  limiting  sensitivity 
of  roughly  one  part  per  thousand  ( 1 016  cm'3)  in  atmospheric  flames.  Under  the  conditions  of 
the  experiments  reported  here,  the  number  density  was  almost  two  orders  of  magnitude  below 
the  value  at  atmospheric  pressure  and  ranged  over  nearly  3.5  orders  of  magnitude  which 
brings  the  lower  limit  to  no  more  than  1014  cm*3.  There  are  several  appropriately  named 
variants  of  the  technique  depending  on  the  geometry  of  the  three  input  beams.  A  review  of 
these  is  given  in  Eckbreth  18  General  reviews  of  the  CARS  technique  can  also  be  found  in 
Arnold  et  al.,19  Tolies,20  Druet21  and  Attal-Tretout.22  The  essential  characteristics  of  the 
technique  are  presented  here.23 


14 


The  CARS  process  is  governed  by  a  fourth-order  susceptibility  tensor,  and  thus  is  a 
member  of  the  nonlinear  family  of  four-wave  mixing  processes.  The  polarization  P(t)  of  a 
scattering  medium  is  related  to  an  incident  electric  field  E(t)  by  the  relation  P(t)  =  e0  x  E(t), 
where  x  is  the  electric  susceptibility  of  the  medium  and  e0  is  the  permittivity  of  free  space. 
Expanding  the  polarization  in  a  power  series  one  obtains 

P{t)  =  (1) 

where  the  x<n>  factors  are  tensors  of  order  n  +  1 .  In  four-wave  mixing,  the  total  applied  field 
is  made  up  of  three  incident  fields  that  can  be  written  as 

£(0  =  ,<V  +£2(co2)e~'“2'  +£,3(co3)e"'“3,+c.c.]  (2) 

where  E/(oJ  is  one  of  the  incident  plane  waves  of  frequency  co;  and  the  propagation 
directions  are  parallel  All  three  terms  in  Eq.  (1)  can,  in  general,  be  simultaneously  nonzero. 
In  a  gas,  however,  the  x0)  term  describes  Rayleigh  scattering  of  the  incident  waves  and  x<2) 
vanishes  identically.  Substituting  Eq.  (2)  into  the  x(3)  term  of  Eq.  (1),  the  general  expression 
for  the  four-wave  mixing  polarization  can  be  written 

^W4=Wm  +  0VWO)=E<>2  2  Xly3^4,0Jm,(0n,Wo)£/((0m)^((0n)£/C0o) 

jkl  (mno) 


(3) 

where  (mno)  indicates  the  sum  over  the  frequency  components  is  restricted  to  the  condition 
indicated  in  the  polarization  argument.  The  number  of  separate  elements  in  x^  is  81 .  The 
number  of  elements  that  are  actually  needed  is  determined  by  symmetry  properties  of  the 
probed  medium.  For  an  isotropic  medium  like  a  gas,  there  are  21  nonzero  elements  that 
possess  the  property  that  any  cartesian  component  that  appears  at  least  once  must  appear  an 
even  number  of  times.  Thus,  because  of  the  equivalence  of  the  three  coordinate  axes,  the  only 
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surviving  elements  are  xxxx,  xxyy,  xyxy,  xyyx  where  xxxx  =  xxyy  +  xyxy  +  xyyx.  In  the 
experiments  reported  here,  the  polarizations  of  the  incident  fields  were  parallel  corresponding 
to  the  xxxx  component. 

The  terms  of  Eq.  (3)  of  interest  here  are  of  the  type  given  in  Eq.  (4),  where  the  first 
expression  corresponds  to  the  case  where  there  are  two  different  pump  waves  at  co,  and  co; 
with  the  probe  wave  at  ©3.  The  second  and  third  expressions  are  essentially  the  same  with 
fo2  =  to,.  (Note  that  the  Vi  factors  are  not  shown  on  the  right  sides  of  these  expressions.) 

P(©4  =  ©j+©2-©3)  =  6x(3)£',£'2£,3* 

P(m4  =  2©,-©3)  =  3  x(3)4e3‘  (4) 

/>(©,  =  2cd2-co3)  =  3  x<3)£22£3‘ 

In  the  case  where  two  incident  frequencies,  one  for  a  pump,  ©p  =  ©,  or  ©2 ,  and  the  other  for 
the  probe,  ©5  =  ©3  (Stokes),  corresponds  to  the  resonance  condition  of  ©v  =  ©p  -  co6 ,  where 
cov  is  the  transition  energy  between  two  vibrational-rotational  levels  of  nitrogen,  then  the 
number  of  terms  contributing  to  the  same  interaction  reduces  to  the  12  terms  summarized  in 
Eq.  (4)  If  the  expression  for  the  electric  field  is  written  with  both  the  temporal  and  spatial 
parts,  thus 


£(©;,r)  =  exp( -/©  /) exp(ik -r)  +  c.c. 


(5) 


and  substituted  in  Eq.  (4)  with  ©2  -  co,  -  cop,  then  the  spatial  phase  of  the  coherent  scattered 
wave,  called  the  anti-Stokes  wave,  at  frequency  co4  =  becomes 

*>  =  (2 kp-ksyr  . 

If  the  incident  waves  are  planar  with  c  and  k  parallel  and  directed  in  the  z  direction,  the 
signal  wave  will  experience  growth  in  this  direction  at  the  frequency  ©„  =  2cop  - 
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Substituting  Eqs.  (4)  and  (5)  into  the  wave  equation,  applying  the  slowly  varying  wave 
approximation  and  defining  the  phase  mismatch  Ak  =  2kp  -  k,  -  one  obtains 


2  ik 


d£(co 


8r 


v  2 

-  =  —ftf3)£2(%)E  ’(<as)exp(i  Akz) 
c 


(6) 


Integrating  Eq.  (6)  over  z  =  0  to  ?,  squaring  the  result  to  obtain  intensity  and  using 


«  ce0 


2 


(7) 


for  each  field,  where  n,  is  the  refractive  index,  Eq  is  the  permittivity  of  free  space  and  c  is  the 
speed  of  light,  the  expression  for  the  C  ARS  signal  beam  intensity  can  be  written 


L  = 


00 


npflsfJas 


c4ej 


A2/2 1  %CARS 1 2  ^SinC  2(Vl  M  0- 


(8) 


where  Xcars  is  the  third  order  susceptibility  for  the  CARS  process.  Maximum  signal  is 
achieved  when  the  phase  matching  condition  Ak  =  0  is  satisfied.  It  should  be  pointed  out  that 
the  pump  intensity  is  on  the  order  of  1  GW/cm2  in  practice  and  the  probe  beam  is  weaker  by 
a  factor  of  roughly  10  but  provides  the  coupling  to  the  resonance  properties  of  the  medium. 
The  strength  of  the  CARS  signal  for  given  peak  intensities  of  the  pump  and  probe  beams 
depends  on  the  square  of  the  CARS  susceptibility,  Xcars  *  which  is  proportional  to  the 

spontaneous  Raman  cross  section  times  the  difference  in  the  population  densities  of  the  two 
levels  connected  by  the  process. 

Diffusion  and  the  Wall  Loss  Flux 

The  probability  of  losing  vibration  quanta  from  nitrogen  molecules  to  a  surface  was 
determined  by  measuring  the  rate  of  decay  of  the  vibrational  state  populations  in  the  center 
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of  the  reactor  tube.  At  a  pressure  of  1 7  Torr  and  a  temperature  of  390  K,  the  V-V  exchange 
rate  is  about  2*  104/s  which  is  roughly  one  tenth  the  rate  molecules  in  the  probe  volume  visit 
the  surface  and  about  103  times  the  wall  relaxation  rate.  There  is,  therefore,  significant 
probability  that  molecules  deactivated  at  the  wall  will  be  re-excited  by  the  V-V  process  before 
returning  to  the  axis  where  the  measurements  are  made.  Thus,  the  observed  rates  of  decay 
of  the  populations  versus  residence  time  (i.e.,  the  length  of  the  specimen  tube  exposed  to  the 
gas  divided  by  the  flow  velocity)  will  be  slower  than  the  actual  wall  rates.  There  are  several 
likely  mechanisms  for  direct  energy  exchange  with  the  surface,  each  possibly  at  a  different 
rate.  Some  mechanisms  result  from  an  event  sequence  referred  to  as  a  channel  wherein  there 
may  be  interactions  between  channels.  There  have  been  numerous  discussions  of  the 
chemistry  and  thermodynamics  on  surfaces  24  25, 26, 27, 28, 29- 30  The  measurements  reported  here 
were  incapable  of  distinguishing  among  these  mechanisms,  but  rather  they  provided  an 
integrated  evaluation  of  all  active  surface  loss  and  gain  mechanisms.  Nevertheless,  the 
accuracy  of  the  measurements  of  these  net  loss  rates  was  not  diminished  by  this  fact. 

Molecules  are  reflected  from  the  wall  without  loss  of  vibrational  energy  with  a 
probability  of  1  -  yv.  With  diffusion,  as  yv  increases,  a  radial  density  profile  in  the  vibrational 
level  populations  develops.  In  the  experiments,  the  length-to-diameter  ratio  of  the  reactor 
tube  was  about  10,  which,  with  the  modest  gas  flow  rates  used,  maintained  a  flat,  “plug-type” 
radial  velocity  profile  along  the  axis  of  the  tube  The  relationship  between  wall  loss  rate  k^^ 

and  probability  yv  is  influenced  by  the  cylindrical  geometry  of  the  reactor  tube.  Because  the 
rate  constant  kM%vaD)  could  be  quite  high,  especially  for  v  =  4  and  5,  the  following  development 
retains  all  corrections  (Note  that  the  “v(wall)”  label  is  dropped  temporarily.)  The  total  flux 
of  particles  (particles  per  unit  area  per  unit  time)  near  a  surface,  J,  is  given  by3' 

J=[N2]j~J{vr)vrdvr  (9) 

where  [N2]  is  the  number  density  of  molecules,  vr  is  the  velocity  normal  to  the  wall  (i.e.,in 
the  radial  direction)  and  f(vr)  is  the  probability  density  of  velocity.  Assuming  the  Maxwell- 
Boltzmann  velocity  distribution  function,  the  flux  becomes 
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(10) 

where  kB  is  Boltzmann's  constant,  T  is  the  temperature,  and  c  is  the  mean  molecular  speed. 
This  expression  must  be  modified  when  there  are  reflection  and  loss  fluxes  at  the  surface.  The 
expression  now  is  given  in  two  equations  which  describe  the  transport  to  and  from  the 
wall.32-”  Thus, 

<u> 

for  vibrational  species  i,  where  the  flux  toward  the  wall  is  J~,  J‘  is  the  flux  in  the  opposing 
direction,  D  is  the  coefficient  of  diffusion  and  w0  is  mean  transport  drift  (flow)  velocity 
towards  the  surface  Since  the  flow  is  only  along  the  tube  axis,  w0  =  0.  Accounting  for  the 
loss  at  the  surface,  we  have  for  that  boundary  condition 

•4™//  =  0  ~Y)4  ~j;  =  0  (12) 

Substituting  Eqs.  (1 1)  in  to  Eq.  (12)  and  simplifying  gives 

D  _  tTO?  r 

*■  4  (i  - 1)  (,3> 

Equation  (13)  shows  that  for  y  =  0  the  flux  reduces  to  that  given  by  Eq.  (10).  In  situations 
where  y  is  small,  the  (1  -  y/2)*1  in  Eq.  (1 3)  can  be  neglected,  but  this  is  not  generally  applicable 
for  the  measurements  reported  here  The  steady  state  diffusion  equation  can  be  written34 

d[N2] 

»0-7 -  =  £V2[tf2]  (14) 

CZ 


j  =  [Ay, 


kBT  j 

T~  = 

2t tm  4 
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where  Uo  is  the  axial  transport  velocity.  A  solution  of  the  form  [N2](r,  z)  =  [N2](0,0)R(r)Z(z) 
is  assumed,  where  (r,z)  are,  respectively,  the  radial  and  axial  coordinates  and  there  is  no 
azimuthal  dependence  Equation  (14)  can  then  be  written  in  the  form  showing  two  equations 
for  R  and  Z,  namely 
S 

_  V8Z  +  D  d2Z  =  _D&R  _  D_dR 
Z  dz  +  Z  dz2  R  dr2  rR  dr  (15) 


ince  the  two  equations  are  independent,  they  must  be  equal  to  a  common  constant  Iq,  where 
the  notation  for  the  species  in  level  v  is  reinstated.  Thus,  we  have 

c*z  _  vaz  _  Kz  _  0 

ar2  D  dz  D 

♦  i*  ♦  .  0  06) 

dr 2  r  dr  D 

The  appropriate  solutions  to  the  Eqs.  (16)  are  Z(z)  =  exp(kzvz)  and  R(r)  =  J0  (Fv  r/r0),  where 
J0  is  the  zero  order  Bessel  function  of  the  first  kind  and  Fv  =  r0  (k,  /D)0'5.  Note  that  kw  and 
lq.  are  wall  rates  in  this  development.  Substituting  the  solution  into  the  Z  equation  and  solving 
for  k^  gives 


+  kJ  in)2 

D\  uQ  ) 


(17) 


In  the  limit  of  lq, «  u02/4D,  Eq.  (17)  reduces  to  kw  =  -lqyu0.  Typically,  this  condition  was 
satisfied  to  better  than  two  orders  of  magnitude.  The  axial  function  now  becomes 

Z(z)  =  exp(  -knz  /  w0)  =  exp(  -kvt)  (18) 
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where  t  =  z/u0  is  the  residence  time  of  the  gas  in  contact  with  the  specimen  wall  Equation 
(18)  describes  the  exponential  decay  of  the  vibrationally-excited  nitrogen  gas  along  the 
specimen  tube  axis  due  to  deactivation  at  the  wall.  Substituting  the  solution  for  the  R 
equation  into  Eq.  (13),  gives  at  the  wall  with  r  =  r0  the  expression 


D—JiW  =  W 


(19) 


where  J,  is  the  first  order  Bessel  function  of  the  first  kind.  If  the  deactivation  probability  yv 
is  known,  then  Eq.  (19)  is  transcendental  in  wall  rate  1^  (i.e.,  k^).  However,  since  the  rate 
was  measured  directly  in  this  work,  Eq.  (19)  can  be  solved  for  yv  giving 

_  4 DFvJx{F)  _  4 DFvJx(Fv) 

'  cr,J0(F,)  (2#> 


In  the  limit  of  small  yv ,  Eq.  (20)  reduces  to  yv  =  2r0kMra„/c.  For  example,  for  typical 
conditions  used  in  this  work,  for  a  value  ofk^,  =  15.5  s'1,  this  approximation  is  8%  smaller 
than  the  value  of  yv  =  5.07x  lO"4  given  by  Eq.  (20). 

Master  Equation  Model 

In  order  to  determine  the  impact  of  the  homogeneous  contributions  to  the  observed 
decay  rates,  a  well-tempered  theoretical  treatment  of  vibrational  energy  exchange  was 
explored  Homogeneous  interactions,  specifically  the  V-V  processes,  mask  the  interactions 
with  the  wall  The  heterogeneous  wall  loss  rates  k^^  determined  in  these  experiments  were 
greater  than  V-T  rates  by  about  three  orders  of  magnitude  while  the  V-V  rates  were  greater 
than  the  wall  rates  by  a  similar  amount.  It  is  commonly  assumed  that  k^,  increases  linearly 
with  v.  Yet,  there  was  some  evidence  early  on  in  this  work  for  deviation  from  a  linear 
relationship  to  a  function  that  peaks  for  v  ^  4.  Moreover,  for  our  conditions,  the  V-V 
“pumping”  process  is  very  efficient  in  creating  a  nonBoltzmann  vibrational  population 
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distribution  in  the  source  gas  that  is  subsequently  altered  by  the  interaction  with  the  wall 
These  issues  indicate  that  there  will  be  a  difference  between  the  observed  decay  of  the 
vibrational  populations  (i.e.,  the  apparent  wall  rate)  and  the  actual  wall  rate.  Clearly,  as 
shown  in  Eq.(20),  the  deactivation  probability,  yv,  is  not  linearly  related  to  due  to  the 
diffusion  process.  Since  the  CARS  probe  is  located  in  the  gas  many  mean  free  paths  from  the 
surface,  it  can  be  expected  also  that  yv  will  not  be  linearly  related  to  the  observed  decay  rate. 
For  these  reasons,  a  concerted  effort  was  made  to  fit  simulations  of  the  decay  of  vibrational 
populations  in  a  cylindrical  tube  based  on  the  Master  Equation  Model  to  all  deactivation  data 
in  order  to  account  for  the  impact  of  homogeneous  interactions  and  diffusion. 

The  simulation  code  was  previously  written  to  predict  the  vibrational  population 
distribution  of  nitrogen  under  only  the  influence  of  loss  terms  from  homogeneous  and 
heterogeneous  vibrational  modes 2  The  code  development  followed  previous  work  of  Bailey35 
and  Capitelli.36  Vibrational  energy  exchange  influences  certain  chemical  reactions  in  plasmas 
and  flames.  In  diatomic  molecules,  this  reactivity  increases  rapidly  with  v.  Theory37  and 
experiment38  show  the  dissociative  attachment  threshold  in  hydrogen  increases  by  five  orders 
of  magnitude  from  v  =  0  to  v  =  6  The  dissociation  of  1%  silane  in  active  nitrogen  was 
measured  to  be  a  factor  of  70  larger  than  expected  from  theory 39  Thus,  a  complete  picture 
of  the  dynamics  of  our  system  requires  more  detailed  knowledge  of  vibrational  pumping  and 
relaxation  processes.  The  V-V  exchange  rates  at  high  v's  are  particularly  difficult  to  measure 
because  the  rates  and  molecular  densities  are  extremely  small.  The  V-V  rates  in  the  lower  v's 
are  better  known  and  are  extrapolated  to  high  v's  in  the  simulation  code  by  the  well  known 
Schwartz,  Slawsky  and  Herzfeld  (SSH)  theory.40 

The  V-V  exchange  rates  are  known  to  be  much  greater  than  the  V-T  rates  at  small  v. 
The  V-T  rate  increases  rapidly  with  v  and  surpasses  V-V  around  v  =  28  depending  on  the 
translational  temperature.  It  is  unlikely  that  the  V-T  interactions  have  any  direct  influence  on 
the  population  distributions  observed  in  this  work,  but  they  were  included  in  the  code  to 
handle  the  upper  population  inversion  extending  to  v  =  47.  The  V-T  process  very  rapidly 
removes  population  for  v  >  47.  Translational  temperatures  in  the  measurements  were  always 
less  than  400  K.  These  temperatures  were  determined  by  recording  a  rotationally-resolved 
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CARS  spectrum  for  the  ground  state  band  (v  =  0  to  v  =  1)  during  each  experiment  and  fitting 
a  theoretical  CARS  spectrum  to  the  measured  spectrum.  The  theoretical  spectrum  was 
generated  by  the  Sandia  CARS  code.41  Rotational  populations  were  observed  to  be 
Boltzmann  distributed  in  all  vibrational  bands,  and  thus  were  equilibrated  to  the  translational 
temperature.  It  was  found  that  the  rotational  temperature  was  independent  of  v  consistent 
with  the  results  of  Massabieaux  et  al.42  It  was  assumed  that  this  held  in  all  experiments 
carried  out  in  this  work.  With  the  modest  temperature  and  low  pressure  used  in  these  studies, 
three  body  collisions  rarely  occur.  Because  the  gas  temperature  was  somewhat  less  than  1000 
K,  the  probabilities  for  multi-quantum  V-V  and  V-T  exchanges  were  very  small  compared  to 
single  quantum  steps,  and  thus  could  be  neglected 43  All  heterogeneous  interactions  also  were 
calculated  on  a  single  quantum  loss  rate  basis.  The  model  assumed  the  molecule  to  be  a 
Morse-type  rigid  oscillator.  These  concepts  and  limitations  were  applied  to  the  Master 
Equation,  which  can  be  written  in  the  generalized  form 


P 


-5* 
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here  is  the  conditional  probability  rate  of  the  transition  m  to  n  and  P„  is  the  occupational 
probability  of  the  property  in  question.44  Here  the  property  is  the  nitrogen  number  density, 
where  each  vibrational  level  is  taken  as  a  separate  species.  The  expression  used  in  the  model 
is  written2 


dNv 

dt 


*,(**,„  +*v,v-l 


(22) 


where  N,  is  the  number  density  (cm'3)  of  a  neutral  nitrogen  molecule  in  the  i*  vibrational  state 
of  the  ground  electronic  state,  the  R^‘s  are  the  exchange  rates  (s’1)  containing  the  V-T  and 
the  V-V  interactions  and  k^^  is  the  wall  deactivation  rate  of  the  Vth  vibrational  state.  The 
rates  in  the  first  two  terms  on  the  right  side  of  Eq.  (22)  can  be  expanded  in  forms  given  by 
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(23) 


*  Ea"  n„ 

k= 0 

P,;.,N  (V, 

*=1 

which  represent  the  V-V  and  V-T  rates  of  single  quanta  energy  exchange  between  molecules 
wherein  population  is  introduced  into  level  v.  For  the  V-V  process,  population  increase  of 
the  v  level  given  in  the  second  term  of  the  first  expression  in  Eq.  (23)  derives  from  one  level 
above  by  transferring  energy  to  a  second  molecule  and  the  corresponding  term  in  the  second 
expression  describes  the  population  increase  from  one  level  below  by  transfer  of  energy  from 
a  second  molecule,  where  the  Qnnpq,  are  the  rate  coefficients  (cm3/s).  The  first  terms  on  the 
right  sides  of  the  two  expressions  in  Eq.(23)  correspond  to  V-T  processes  which  increase  the 
population  of  the  v  level,  where  the  Pm  „  are  the  rate  coefficients  (cm3/s)  and  N  is  the  total 
number  density  of  the  gas.  The  corresponding  loss  rates  given  by  the  third  and  fourth  terms 
of  Eq.(22)  have  terms  analogous  to  Eq.(23). 
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ANALYSIS 

Data  Collection.  Corrections  and  Normalization 

Plots  of  typical  spectra  for  the  v  =  0,  1,2  and  3  bands  are  shown  in  Fig.  3  The 
measurements  were  made  at  the  peaks  of  the  J  =  8  lines  of  the  vibrational  bands.  This 
simplified  both  the  deactivation  analyses  and  the  saturation  corrections;  however,  because  of 
the  narrow  width  of  the  rotational  lines,  the  CARS  signals  were  extremely  sensitive  to  the 
wavelength  setting  of  the  dye  laser.  Since  the  YAG  laser  was  operating  single  mode  (0.006 
cm'1),  the  observed  width  was  the  convolution  of  the  Doppler  lineshape  of  the  rotational  line 
with  the  lineshape  of  the  dye  laser.  Since  the  Doppler  width  was  0.006  cm*1  and  the  dye  laser 
linewidth  was  about  0.1  cm'1,  the  observed  width  was  determined  by  the  dye  laser.  The 
smallest  wavelength  step  provided  by  the  dye  laser  was  0.001  nm  (0.027  cm'1),  however,  the 
resetability  was  two  to  four  times  larger  than  the  smallest  step.  Thus,  it  was  very  difficult  to 
dependably  move  from  one  vibrational  band  to  the  next  without  misalignment  of  the  dye  laser 
wavelength.  Therefore,  the  rotational  line  in  each  vibrational  bands  was  scanned  and  then 
fitted  to  the  squared  Lorentzian  function.  This  procedure  provided  a  precise  measurement 
of  the  peak  CARS  intensity  with  a  much  relaxed  requirement  for  the  initial  wavelength  setting 
of  the  dye  laser  This  type  spectral  measurement  is  termed  “skip  spectra.” 

Skip  spectra  were  obtained  in  both  deactivation  and  saturation  experiments.  In  a 
typical  experiment,  which  involves  several  runs  over  the  desired  vibrational  bands,  each  J  = 
8  line  is  scanned  with  16  spectral  points  At  each  spectral  point  in  a  scan,  the  average  of  200 
to  500  laser  pulses  was  recorded.  For  each  laser  pulse,  the  PMT  produces  a  signal  CARS 
pulse  that  is  integrated  with  1 0-bit  precision  by  the  LeCroy  ADC.  The  signals  from  the  two 
laser-monitoring  PIN  diodes  are  also  integrated  with  the  ADC  and,  with  the  integrated  CARS 
signals,  are  averaged  at  each  spectral  point.  Signal  data  are  recorded  in  double  precision  for 
each  spectral  point  and  stored  in  files  as  sums.  This  includes  the  sum  and  the  sum  of  squares 
for  the  CARS,  dye  and  YAG  PEN  diode  signal  as  well  as  the  wavelength  and  number  of  pulses 
summed.  For  the  laser  covariance,  the  sums  of  sums  and  sums  of  differences  of  diode  signals 
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CARS  Signal  (a.u.)  CARS  Signal  (a.u.) 


Raman  Shift  (cm1) 


Figure  3.  Plots  of  the  v  =  0,  1,  2  and  3  rotational- vibrational  Q  branch  bands 
of  N2  obtained  with  the  CARS  system  in  Fig.  2. 
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and  their  sum  of  squares  are  also  written  to  the  file.  Separately,  hand  records  were  kept  of 
the  optical  attenuation  of  the  CARS  beam,  electronic  channel  gains  and  the  beam  powers 
between  each  run.  The  files  and  notes  were  used  to  determine  the  rescaled  spectral  value, 
saturation  correction  and  the  statistics 

Rescaling  of  data  differs  between  a  saturation  and  a  deactivation  experiment.  For  the 
saturation  case,  each  band  v  is  normalized  separately  to  known  YAG  and  dye  laser  powers. 
In  the  deactivation  case,  an  entire  run  is  rescaled  using  the  average  PIN  diode  signals  from 
the  v  =  0  measurements  as  the  reference  values  and  then  the  raw  CARS  signals  are  corrected 
for  saturation.  The  rescaling  places  the  measurements  on  each  vibrational  band  in  a  given 
deactivation  experiment  on  the  same  the  basis  of  YAG  and  dye  powers  so  that  the  relative 
band  populations  can  be  extracted.  Examples  of  these  calculations  are 
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where  the  sums  are  for  each  spectral  point  i  of  the  band  v,  Scj  is  the  CARS  signal  sum,  S2ci 
is  the  sum  of  the  squares  of  the  CARS  signals,  Yf  is  the  sum  of  the  yag  PIN  diode  signals,  and 
Dj  is  the  sum  of  the  dye  PIN  diode  signals  The  lower  case  quantities  are  the  pulse  values  and 
the  means  over  npts  pulses  are  calculated  later 

Each  rotational  line  was  scanned  with  np  points  that  includes  two  points  in  the 
background.  The  CARS  background  signals  were  measured  close  to  the  band  heads  of  each 
vibrational  band  in  regions  essentially  free  of  any  detectable  rotational  lines.  The  background 
measurement  is  subtracted  from  the  average  signals,  which  removes  any  electronic  offset  and 
nonresonant  CARS  signal  The  CARS  channel  offset  came  from  two  external  sources.  One 
was  a  slight  sensitivity  to  extraneous  electrical  emissions,  principally  from  the  flash  lamp 
power  supply,  and  the  other  was  leakage  of  the  532  nm  light  onto  the  PMT.  Electronic 
offsets  in  PIN  diode  channels  did  not  change  more  than  five  counts  in  three  years  and  was, 
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thus  constant  throughout  an  experiment.  Offsets  in  the  PIN  diode  channels,  Doff  and  Yoff , 
were  determined  by  observing  the  dark  count  value  for  each  detector,  which  was  typically  25 
counts. 

CARS  background  exhibited  a  slight  gradual  change  across  a  vibrational  band.  This 
characteristic  was  constant  through  an  experiment  and  probably  over  many  experiments. 
Actually,  the  nonresonant  background  was  weak  or  nonexistent.  A  conventional  background 
correction  uses  measurements  on  either  side  of  the  feature  to  predict  the  background  at  the 
peak  of  the  line.  In  this  work,  the  value  on  one  side  was  subtracted,  which  provided  a 
consistent  and  expedient  method.  The  error  was  usually  well  under  10  counts  or  1%.  This 
error  was  largest  in  the  high  hot  bands  where  the  signal-to-background  ratio  was  the  lowest. 


The  band  averages  for  the  dye  laser  and  YAG  laser  PIN  diode  signals  over  np  spectral 
points  are  computed  from 
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where  Doff  and  Yoff  were  25,  and  the  rescaled  CARS  signal,  Sc2i ,  was  calculated  for  each 
spectral  point  i  in  band  v  as  given  in  by 
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which  normalizes  each  spectral  datum  to  the  YAG  and  dye  power  averages  for  that  scan.  The 
variances  of  the  means  of  the  CARS  signals,  the  dye  and  YAG  PIN  diode  signals  and  the 
covariance  between  the  dye  and  the  YAG  diode  signals  were  also  calculated.  The  correlation 
coefficient  between  diode  signals  was  typically  between  0.2  and  0. 1  or  less  which  supports 
the  normalization  using  the  averages 

Saturation  corrections,  derived  from  calibration  experiments,15  were  used  to  correct 
the  reduction  in  response  caused  by  stimulated  Raman  scattering  (SRS)  which  pumps 
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population  directly  into  the  upper  level  simultaneously  with  the  CARS  process.15-  45 
Saturation  calibration  measurements,  which  gave  the  relationship  between  the  peak  CARS 
signals  and  laser  powers,  would  be  conducted  typically  within  no  more  than  20  days  of  the 
deactivation  experiments  to  which  the  calibration  was  applied  providing  that  the  optical 
alignment  of  the  system  had  not  been  disturbed  The  average  dye  laser  power  was  varied  in 
steps  from  ~33  to  ~  10  mW  for  the  saturation  measurements.  Usually,  30  mW  (1.5  mJ/pulse 
at  20  Hz)  was  the  power  used  in  the  deactivation  experiments.  Power  incident  on  the  main 
chamber  was  a  convenient  unit  for  the  calibration  of  two  PIN  diodes.  The  J  =  8  line  was 
scanned  in  at  least  three  of  the  normally  detected  hot  bands  of  v  =  1,  2,  3,  and  4  at  each  dye 
laser  power  step  while  holding  all  other  system  variables  constant.  Nine  to  12  power  steps 
were  recorded.  Using  the  procedure  described  by  Yaney  and  Parish,15  the  average  CARS 
signal  (PMT  counts)  in  band  v,  Sc2v,  at  the  J  =  8  peak  and  the  laser  powers  Yv  and  Dv  (PIN 
diode  counts)  were  combined  to  plot  (Se2v/Dv  -Yv2)05  versus  Dv  Yv  for  each  v-band. 
Parameterized  in  this  way,  the  data  gives  a  linear  relationship  for  each  vibrational  band  with 
an  intercept  bv  at  zero  laser  power  and  a  slope  n^.  as  given  in 
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where  n\  is  negative  The  square  root  of  the  ratio  of  the  unsaturated  to  saturated  signal,  Ssat  v 
gives  the  correction  factor  for  each  vibrational  level  as15 
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The  ratio  myb,,  in  Eq.(28)  generally  is  a  linear  function  of  v.15  Therefore,  doing  a  linear 
regression  of  these  ratios  to  v  +  1  from  as  few  as  three  v  bands  provided  a  convenient  way 
to  accomplish  the  correction  of  the  deactivation  data  for  all  v  bands  under  each  particular  set 
of  experimental  circumstances.  Since  the  PIN  diode  responses  to  the  powers  of  their 
respective  lasers  had  been  calibrated,  the  saturation  calibrations  could  be  used  to  correct  the 


29 


deactivation  data  using  count  values  in  Eq.(28).  It  can  be  shown  that  the  methods  used  are 
equivalent  to  using  power  values.  The  deactivation  and  saturation  experiments  were 
connected  by  diode  counts-per-mW  (cnts/mW)  ratios  which  were  used  to  convert 
deactivation  power  readings  to  appropriate  saturation  values. 

The  saturation  measurements  were  normally  carried  out  by  varying  the  dye  laser 
power  by  rotating  a  halfwave  plate  just  in  front  of  a  polarizer  while  holding  the  YAG  power 
constant.  To  check  that  the  saturation  effect  depended  on  the  product  of  the  two  laser 
intensities  in  the  probe  volume  as  shown  in  Eq.(28),  the  YAG  power  was  varied  in  one 
experiment  while  holding  the  dye  power  constant.  The  CARS  signal  had  the  same  saturation 
dependence  on  the  incident  YAG  power  as  it  did  by  varying  the  dye  power. 

The  data  at  the  spectral  points,  Se2,  are  free  of  offset,  background  and  any  laser  power 
drift  during  the  scan  over  the  line.  These  spectral  data  of  the  rotational  lines  were  fitted  to 
a  Lorentzian  squared  function  to  obtain  the  peak  heights  as  discussed  in  the  next  section.  The 
rotational  peak  heights  were  then  corrected  for  saturation  and  the  resulting  values  were 
rescaled  to  the  average  YAG  and  dye  laser  powers  present  during  the  measurements  on  the 
v  =  0  to  1  ground-state  band  in  that  experiment.  However,  the  optical  attenuation  and  the 
amplifier  gain  were  adjusted  for  each  hot  band  to  fit  the  observed  signals  well  within  the 
1 024-count  digitization  range  of  the  ADC.  The  dynamic  range  of  signals  over  the  four  or  five 
vibrational  bands  was  greater  than  7  *  105.  This  dynamic  range  was  covered  by  the  neutral 
density  filter  wheel  that  provides  attenuations  from  unity  to  2145  and  the  amplifier  that  has 
10  gain  steps  from  unity  to  10.  Generally,  the  attenuation  and  the  gain  were  adjusted  to  keep 
the  maximum  CARS  signal  pulses  within  900  counts  of  the  ADC.  The  CARS  signals 
displayed  considerable  fluctuations  attributed  to  the  mode  beating  in  the  dye  laser  pulse  which 
made  these  adjustments  very  critical  to  acquiring  viable  data. 

Fitting  the  rotational  line 

In  order  to  achieve  a  high  degree  of  precision  in  determining  the  CARS  signal 
strength,  the  J  =  8  rotational  line  in  each  vibrational  band  was  measured  by  scanning  this  line 
using,  typically,  14  steps  of  0.001  nm  (0.027  cm'1)  and  200  laser  pulses  per  step  plus  two 
background  steps  The  background  signal  level  was  removed  from  these  raw  data  and  the 
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square  of  the  Lorentzian  profile  was  used  to  fit  to  the  background-free  data  using  nonlinear 
regression  analysis.  This  profile  can  be  written  for  narrow  lines  in  wavelength  form  as 


L(k)  = 


w\ 


l+[(w2-X,)/w3]2 


(29) 


where  w,  is  the  square  root  of  the  peak  CARS  amplitude  of  the  line,  w2  is  the  wavelength  at 
the  peak,  w3  is  the  linewidth  (FWHM)  and  A,,  is  the  wavelength  at  the  i*  point  of  the  scan. 
Because  any  background  signal  was  removed  prior  to  applying  Eq.  (29),  no  additive  constant 
was  needed  in  the  expression  Parameter  Wj  is  proportional  to  the  population  difference 
N(v,  J  =  8)  -  N(v  +  1,  J  =  8),  which  was  then  used  to  compute  the  relative  populations  among 
the  observed  vibrational  bands.  The  fitted  values  for  w,  were  recorded  in  a  disk  file  together 
with  the  PIN  diode  count  averages,  Dv  and  Yv .  These  values  of  w,  were  then  fully  corrected, 
as  stated  at  the  end  of  the  last  section,  to  the  ground  state  laser  powers  and  for  saturation. 
The  corrected  values,  wp.,  for  band  v  are  obtained  from 
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here  Y0  and  D0  correspond  to  the  average  ground  state  (v  =  0  to  v  =  1)  laser  power  values, 
G  is  the  amplifier  gain  and  A  is  the  optical  attenuation  used  in  the  band  v  measurements. 

Obtaining  the  observed  deactivation  rates 

Following  the  various  normalizations  and  corrections  of  the  peak  values  of  the  J  =  8 
line  in  each  vibrational  band,  the  data  for  the  v  bands  at  each  residence  time  were  processed 
together  to  extract  fractional  populations.  The  residence  time  is  simply  the  length  of  the 
specimen  tube  exposed  to  the  vibrationally-excited  gas  divided  by  the  gas  flow  velocity  in  the 
tube  It  is  important  to  emphasize  that  the  peak  line  strengths,  w^  ,  obtained  at  a  given 
residence  time  were  processed  together  independent  of  the  data  at  other  residence  times.  As 
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shown  in  Eq.(8),  the  line  strength  is  proportional  to  the  CARS  susceptibility,  Xcars,  which  is 
given  by18 
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where  the  frequency  <%  is  the  frequency  of  the  transition  between  levels  a  and  b,  do/dQ  is  the 
differential  cross  section  of  the  linear  Raman  process  which  depends  on  \,p  and  5  identify 
pump  and  probe  (Stokes)  quantities,  h  is  Planck’s  constant  divided  by  2ji,  Na  and  Nb  are  the 
population  densities  and  Tab  is  the  line  width.  The  total  susceptibility  includes  the  real  and 
imaginary  parts  obtained  from  Eq.(3 1 )  and  a  nonresonant  component  as  shown  in  Eq.(32). 

Xtor  “  'Xcars  +  ‘Xcars  +  X^r  (32) 

If  the  nonresonant  component  has  sufficient  magnitude  it  can  complicate  the  functional  form 
of  the  CARS  spectrum  due  to  the  squaring  in  Eq.(8).  In  this  work,  was  essentially  zero 
At  resonance  when  cop  -  ws  =  co^,  which  corresponds  to  the  peak  of  the  line,  only  the 
imaginary  part  of  Eq.(3 1 )  remains.  Thus,  using  the  peak  line  strength  simplifies  the  analysis. 
It  is  assumed  in  this  analysis  that  the  rotational  temperature,  and  therefore  the  rotational 
population  distribution,  is  equilibrated  with  the  translational  or  gas  temperature.  The 
population  distribution  over  v  and  J  is  adequately  described  by 
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where  N  is  the  total  density  of  molecular  nitrogen,  Z  is  the  partition  function,  G(v)  is  the 
vibrational  term  energy  from  the  ground  state,  v  and  J  are,  respectively,  the  vibrational  and 
rotational  quantum  numbers,  Be  is  the  equilibrium  rotational  constant,  ae  is  the  vibration- 
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rotation  interaction  constant,  Tv  is  the  vibrational  temperature  of  band  v  corresponding  to  the 
v  to  v  +  1  transition,  Tr  is  the  rotational  temperature,  k  is  the  Boltzmann  constant  and  g,  is  the 
statistical  weight  due  to  nuclear  spin  degeneracy  which  depends  on  J  for  nonzero  nuclear  spin. 
The  partition  function  is  found  by  summing  Eq.(33)  over  v  and  J;  however,  it  can  be  factored 
into  vibrational  and  rotational  parts  with,  at  most,  a  1%  error  by  neglecting  the  right  hand 
factor  involving  v  and  J  in  Eq.(33).  Thus, 
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and  Z,  is  the  sum  of  the  J-only-dependent  factors  in  Eq.(33).  The  population  difference  for 
a  Q  Branch  transition  where  AJ  =  0,  ANViJ  =  NV>J  -  Nv+W  can  be  written  as 

WvJ  =  K^NvKvJ  -  Nv+lKvtU) 
where  KvJ  =  exp[ae(v  +  0.5)J(J+l)/^r  ]  (3g) 

and  K'j  =  gj(2J+l)exp[-BeJ(J+\)/kTr]/Zr 


From  Eqs.(8)  and  (33),  the  square  root  of  the  CARS  signal  can  be  written 
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where  cv  is  the  Raman  cross  section  of  the  v  vibrational  band  and  h(v)  corrects  for  the 
deviation  of  the  Morse  potential  from  an  harmonic  oscillator  .46  The  width  of  the  Doppler- 
broadened  rotational-vibrational  line  depends  on  the  frequency  of  the  transition  and  the  gas 
temperature.  This  introduces  another  correction  factor,  Dv ,  given  by46 
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The  strength  of  the  transition  is  inversely  proportional  to  the  width.  E^. ,  is  the  energy  of  the 
rotational-vibrational  level  and  Eq.(37)  references  the  strength  of  a  transition  in  band  v  to  the 
ground  state  band  v  =  0.  All  measurements  were  made  on  J  =  8.  Combining  Eqs.(35)-(37) 
assuming  constant  laser  powers,  perfect  phase  matching  and  constant  beam  intersection 
geometry,  the  square  root  of  the  CARS  signal  becomes 


Sv.j  =  Cj  (v  H- 1)  hvDv(KvJNv  -  K^u  NVJ 


(38) 


where  C,  contains  the  laser  power  factors  and  all  the  constants  in  Eqs.  (8),  (31),  and  (34-(36). 
At  some  v  band  the  population  will  be  so  small  that  the  signal  will  be  unmeasurable.  It  is 
assumed  that  populations  at  higher  v’s  had  negligible  affect  on  the  measured  populations  in 
the  lower  v  bands.  This  permitted  a  truncation  of  the  calculation  for  extraction  of  fractional 
populations.  Thus,  by  solving  Eq.(38)  for  Nv  gives  a  summation  from  v  to  the  maximum  v 
measured,  v^ ,  which  for  J  =  8,  yields  Eq.(39),  where  Nv’  =  NvCj  is  a  relative  population. 
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In  this  summation,  the  population  of  the  vmax  +1  band  is  set  to  zero.  The  total  relative 
population,  N’,  is  then  the  sum  of  Eq.(39)  from  v  =  0  to  v  =  v^.  Thus,  the  fractional 
populations,  Fv,  are  found  from 
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Various  schemes  were  explored  to  extrapolate  the  calculation  of  Fv  to  include  higher  values 
of  vmax  based  on  theoretical  projections.  Although  these  extrapolations  produced  small 
changes  in  the  values  of  F4  and  Fs  compared  to  the  above  truncation  procedure,  these  changes 
had,  at  most,  only  a  small  impact  (<10%)  on  the  observed  relaxation  rates. 
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DATA  AND  RESULTS 


Tables  1-6  summarize  the  results  of  this  project.  Reproducibility  of  these  results 
varied  because  of  accumulated  surface  changes  due  to  use  of  the  specimen  in  many 
experiments,  planned  treatments  of  the  specimens,  undetermined  variations  in  the  specimens 
or  uncontrollable  variations  in  the  experimental  system.  Several  experiments  with  Pyrex 
tested  these  possibilities.  Tables  1,  3  and  5  report  the  observed  relaxation  rates,  k^,  and 
Tables  2,  4  and  6  give  the  corresponding  apparent  deactivation  coefficients,  yw  computed 
directly  from  the  k,.  (viz.,  instead  of  k^)  values  using  Eq.(20).  Generally,  k^  is  expected 
to  increase  for  higher  pressure  and  decrease  with  increased  temperature.  For  a  given 
pressure,  the  temperature  for  specified  discharge  conditions  was  determined  by  the  gas  flow 
rate  Higher  flow  rates  gave  higher  temperatures  due  to  increased  heat  transport.  The 
observed  relaxation  rate,  however,  is  a  combination  of  the  interactions  of  the  molecules  with 
the  wall  and  the  interactions  with  other  molecules  through  the  V-V  exchange  process.  This 
greatly  complicates  the  interpretation  of  the  data.  The  fitting  of  the  Master  Equation  Model 
attempts  to  extract  the  wall  relaxation  rate,  but  the  interplay  between  these  two  interactions 
can  vary  significantly  from  experiment  to  experiment.  The  results  of  these  fittings  are 
discussed  in  the  next  section. 

Except  for  gold,  all  surfaces  were  found  to  be  nitrided  by  the  afterglow  due  to  the 
small  density  of  nitrogen  atoms.47  This  density  was  estimated  to  be  0.55%  of  [N2]  from 
measurements  of  the  decay  of  the  first  positive  emission,  B3FIg  to  A3  u,  versus  residence  time. 
Surface  measurements  on  Pyrex  using  XPS  showed  nitrogen  bonded  principally  to  the  silicon 
and  secondarily  to  the  oxygen  surface  atoms  of  witness  coupons.48  The  complete  XPS  results 
are  reported  in  a  later  section.  Coverage  was  determined  to  be  between  8%  and  10%  of  the 
surface  with  one  case  showing  about  30%.  It  is  clear  that  exposure  to  the  afterglow  changes 
the  surface  in  some  way.  The  details  of  the  bonding  mechanisms  cannot  be  obtained  nor 
completely  understood  from  these  experiments.  Nitrogen  adatoms  are  expected  to  make  the 
surface  passive  by  permanently  occupying  surface  adsorption  wells  generated  by  dangling 
bonds  They  also  change  the  principal  surface  resonance  vibration  away  from  that  of  the 
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virgin  surface  which  could  provide  more  efficient  coupling  to  vibrating  molecules.  Thus, 
adatoms  may  either  increase  or  reduce  the  surface  deactivation  rate  and,  in  the  same  respect, 
the  number  of  surface  collisions  a  molecule  makes  or  the  number  of  quanta  lost  per  collision 
may  be  impacted  by  the  nitriding. 


Pvrex 

Table  1  gives  the  observed  values  ofk,  with  fitting  uncertainties  and  Table  2  gives  the 
corresponding  apparent  deactivation  coefficients  y^.  These  tables  also  contain  results  of 
Master  Equation  modeling  efforts.  In  Table  1,  “nitrided”  indicates  the  surface  had  been 
exposed  to  active  nitrogen  (i.e.,  afterglow)  in  several  previous  experiments  and  was  found 
from  XPS  measurements  to  have  nitrogen  bonded  to  the  surface;  “neat”  refers  to  the  first  time 
this  specimen  was  used,  taken  as  is  from  the  glass  shop  and  treated  only  to  a  general  cleaning; 
“+column”  refers  to  an  experiment  where  the  tube  was  placed  in  the  positive  column  for  33 .7 
hours  before  the  deactivation  measurements;  and  “exposed”  refers  to  an  experiment  where 
the  tube  spent  39.5  hours  in  the  afterglow  before  deactivation  measurements.  The  cleaning 
procedure  generally  consisted  of  ultrasonic  degreasing  in  methylene  chloride  followed  by  ion 
cleaning  at  about  5  Torr  in  flowing  argon  for  30  to  60  minutes. 

Experiments  1  and  2  used  the  same  tube  wherein  the  latter  experiment  occurred  after 
additional  exposure  of  the  tube  to  the  afterglow  thereby  increasing  the  degree  of  nitriding 
It  would  seem  that  the  additional  nitriding  in  Experiment  2  is  responsible  for  the  higher  values 
°f  Yva-  The  lower  flow  rate  and  the  resulting  lower  temperature  introduces  ambiguity  in  any 
attempt  to  interpret  these  data.  Changes  due  to  exposure  to  the  afterglow  were  also  reported 
by  Black  et  al  using  a  microwave  source8;  however,  they  report  a  decrease  of  y,  for  Pyrex 
after  a  24  hour  exposure  to  the  afterglow.  They  found  no  change  with  exposure  when  using 
a  thermal  source.  This  can  be  explained  by  the  fact  that  they  measured  only  the  v  =  1  band 
using  a  1-m  long  specimen  tube.  With  the  microwave  source,  as  with  the  positive  column 
source  used  in  this  work,  the  vibrational  population  is  distributed  across  all  low-v  states 
corresponding  to  a  vibrational  temperature  typically  well  above  the  2000  K  they  report  for 
their  thermal  source.  Moreover,  there  is  no  surface  nitriding  using  a  thermal  source.  Thus, 
measuring  only  the  v  =  1  population  produced  by  a  discharge  source  is  not  sufficient  to 
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correctly  characterize  the  population  changes  due  to  wall  relaxation  over  long  residence 
times,  especially  with  V-V  pumping  and  nitriding  effects. 

Experiments  3  and  4  were  carried  out  to  check  for  reproducibility  and  the  degree  of 
dependence  on  pressure.  The  data  and  fitting  results  for  Experiment  3  are  given  in  Fig.  4. 
The  flow  rate  in  Experiment  4  was  not  increased  to  compensate  for  the  higher  pressure; 
therefore,  the  temperature  in  Experiment  4  was  not  the  same  as  in  Experiment  3. 
Nevertheless,  the  relaxation  rates  in  these  two  experiments  are  the  same  within  fitting 
uncertainties  Experiment  5  attempted  to  achieve  a  high  degree  of  nitriding.  Although  the 
Yva  values  are  larger  than  the  values  obtained  for  the  “neat”  specimens,  the  rate  of  increase  of 
Yva  with  v  for  the  “+  column”  specimen  is  about  half  that  of  the  neat  specimens.  Experiments 
4,  6  and  7  used  one  tube  and  Experiment  7  followed  the  next  day  after  Experiment  6 
Experiments  6  and  7  were  run  to  test  reproducibility.  The  yva  values  show  reasonable 
agreement  within  the  fitting  uncertainties  plus  the  estimated  5  to  10%  uncertainty  in  the 
repeatability  of  the  set  up  of  the  experimental  conditions.  Comparison  of  Experiments  6  and 
7  to  Experiment  4  again  shows  the  rate  of  increase  of  with  v  to  be  slower  after  the  tube 
was  exposed  to  the  afterglow  An  experiment  was  attempted  at  a  pressure  of  0.5  Torr  also 
on  the  this  specimen  Although  the  signals  were  very  weak,  the  observed  relaxation  rates  for 
v  =  1  and  2  were  about  half  the  values  obtained  in  Experiments  6  and  7  confirming  the 
expected  pressure  dependence. 


Gold  Foil 

The  results  reported  in  Tables  2  and  3  are  for  the  same  gold  foil  and  are  tabulated  in 
order  of  acquisition.  The  foil  was  formed  to  completely  cover  the  inside  of  a  Pyrex  tube. 
Gold  was  chosen  because  surface  changes  due  to  reaction  with  active  nitrogen  were  not 
expected.  Thus,  experiments  could  be  conducted  to  assess  the  effects  of  temperature 
(determined  by  the  gas  flow  rate)  and  pressure  in  the  absence  of  nitriding.  It  should  be  noted, 
however,  that  the  surface  did  change  from  a  calendered  gloss  to  a  slightly  darker,  very  fine 
buff  This  surface  roughening  could  only  be  caused  by  collisions  of  vibrationally  hot  nitrogen 
during  the  experiments  or  by  the  argon  positive  column  used  for  surface  cleaning.  The  first 
three  experiments  in  Tables  3  and  4  used  three  different  gas  flow  rates  in  order  to  vary  the  gas 
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Figure  4.  Relaxation  measurements  on  neat  Pyrex  at  360  K  and  1 1  Torr:  the  points  are  the  data; 
the  solid  lines  are  exponential  fits  to  the  data  with  the  values  given;  the  dashed  lines  are  the 
results  of  the  fitting  of  the  Master  Equation  Model  to  the  data;  and  the  vibrational  temperature 
of  the  v  =  1  band,  T, ,  is  given  by  the  triangles  joined  by  the  cubic  spline  curve. 


Vibrational  Temperature  T1  (K) 


temperature.  The  tube  was  mounted  in  the  chamber  using  glass  beads  that  provided  good 
thermal  isolation.  Thus,  it  is  assumed  that  the  tube  temperature  was  equal  to  the  gas 
temperature  as  measured  from  a  CARS  rotational  spectrum.  As  expected,  the  relaxation  rates 
decreased  uniformly  with  increase  in  temperature.  The  last  three  experiments  were  intended 
for  pressure  studies  Comparing  Experiments  1  and  4,  which  are  at  the  same  temperature, 
show  significant  decreases  in  relaxation  rates  and  the  corresponding  y„,  values  with  drop  in 
pressure;  however,  the  opposite  change  is  shown  in  Experiments  2  and  5  which  were 
performed  at  a  higher  temperature.  It  is  likely  that  this  apparent  contradiction  is  a 
consequence  of  the  complexities  of  the  mixing  of  the  wall  relaxation  effects  and  the  V-V 
interactions.  The  data  and  fitting  results  for  Experiment  2  are  given  in  Fig.  5. 

Selected  Materials 

Tables  5  and  6  report  the  results  on  a  variety  of  materials  in  the  same  form  as  the 
previous  tables.  The  two  stainless  steel  specimens  showed  the  highest  relaxation  rates  as 
would  be  expected  and  as  reported  by  Black  et  al.8  The  304  alloy  was  measured  at  the  start 
of  these  studies  when  the  discharge  tube  was  operated  at  a  higher  temperature;  however,  the 
inside  surface  of  this  tube  had  been  machined  to  a  fairly  rough  surface.  On  the  other  hand, 
the  AMS  3 1 2  tube  had  a  very  smooth  “drawn”  surface.  Although  the  alloys  are  different,  the 
high  relaxation  rates  obtained  with  the  304  specimen  is  attributed  to  the  large  effective  surface 
area  Figure  6  shows  plots  of  the  data  and  the  fitting  results  for  AMS  312.  Two  sets  of 
values  are  given  for  aluminum  which  represent  experiments  on  two  different  tubes  of  the  same 
alloy.  The  results  are  in  reasonable  agreement  which,  again,  shows  that  the  measurement 
procedures  were  carried  out  in  a  consistent  manner.  The  results  for  titanium  were  a  surprise. 
The  same  tube  was  used  in  both  experiments.  Experiment  3  was  first  and  showed  virtually 
no  relaxation  of  the  v  =  1  band  while  Experiment  4  showed  that  the  surface  had  been  altered 
by  the  exposure  to  active  nitrogen.  This  nitriding  again  reduced  the  rate  of  increase  of 
relaxation  rate  with  v,  which  is  a  factor  of  two  between  Experiments  3  and  4.  Titanium  gave 
the  lowest  relaxation  rates  observed  in  these  studies,  which  has  to  be  attributed  to  the  oxide 
layer  on  the  titanium  surface.  Teflon,  chosen  to  represent  polymer-based  materials,  was 
expected  to  be  benign  This  material,  indeed,  gave  the  lowest  relaxation  rates 
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Fractional  Population 


Figure  5.  Relaxation  measurements  on  gold  foil  at  351  K  and  17  Tort:  the  points  are  the  data; 
the  solid  lines  are  exponential  fits  to  the  data  with  the  values  given;  the  dashed  lines  are  the 
results  of  the  fitting  of  the  Master  Equation  Model  to  the  data;  and  the  vibrational  temperature 
of  the  v  =  1  band,  T,  ,  is  given  by  the  triangles  joined  by  the  cubic  spline  curve. 
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Figure  6.  Relaxation  measurements  on  AMS  312  stainless  steel  at  340  K  and  14  Torr:  the  points 
are  the  data,  the  solid  lines  are  exponential  fits  to  the  data  with  the  values  given;  the  dashed 
lines  are  the  results  of  the  fitting  of  the  Master  Equation  Model  to  the  data;  and  the  vibrational 
temperature  of  the  v  =  1  band,  T, ,  is  given  by  the  triangles  joined  by  the  cubic  spline  curve. 
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except  for  titanium.  With  two  fluorine  atoms  per  carbon  atom.  Teflon  contains  a  high  fluorine 
fraction  which  would  otherwise  contribute  to  strong  relaxation.  The  stability  of  these  bonds 
is  generally  quite  high,  and  therefore  Teflon  would  be  expected  to  have  only  a  weak 
interaction  with  vibrational-excited  species.  Nevertheless,  it  was  found  very  important  to 
keep  the  gas  temperature  low  with  this  material,  otherwise  the  surface  would  be  altered  and 
release  contaminates  into  the  flow. 


Table  1 :  Observed  relaxation  rates  k..  for  Pyrex 


The  %uncertainty  in  yva  *  (1+0.01  kj  x%uncertainty  in 


Table  5.  Observed  relaxation  rates  k,.  of  selected  materials 


Exp 

Material 

H 

■ 

■ 

1 

H 

k4 

(s'1) 

1 

Stainless  steel 

460 

17 

1000 

53.1 

108 

138 

S  304 

II 

2 

Stainless  steel 

340 

14 

800 

19 

44 

70 

135 

AMS  312 

±9% 

±5% 

±6% 

±7% 

3 

Aluminum 

392 

17 

1000 

12 

o 

31 

38 

3003 

349 

14 

800 

8.3 

El 

32 

42 

HI 

±7% 

±7% 

±6% 

±6% 

D 

Titanium 

382 

14 

1100 

«0 

5 

8 

13 

■ 

AMS  4943 D 

±26% 

±21% 

±20% 

5 

Titanium 

375 

14 

800 

1.3 

ESI 

6.4 

8.1 

AMS  4943D 

±22% 

Wm 

±9% 

±8% 

6 

Teflon 

334 

14 

800 

EBI 

8.8 

14 

20 

IHyJ 

±10% 

±7% 

±8% 

TabIe  6  Apparent  deactivation  coefficients  YVJ  x  104  of  selected  materials' 


Exp 

Surface 

Yu 

Y2a 

Y* 

Y4a 

1 

SS  304 

18 

53 

96 

— 

2 

SS  312 

7.8 

23 

58 

104 

3 

A1  3003 

4 

7.0 

12 

15 

2.8 

7.3 

12.7 

18.3 

4 

Ti  4943D 

*0 

1.5 

2.5 

4.3 

5 

Ti  4943 D 

0.4 

1.2 

2.1 

2.9 

6 

Teflon 

1.1 

2.8 

4.7 

7 

The  %uncertainty  in  *  (l+O.Olk^/ouncertainty  in  k. 


Master  Equation  Modeling 

The  values  given  in  Tables  2,  4  and  6  include  the  influence  of  the  V-V  processes  as  well  as 
interactions  with  the  wall.  The  Master  Equation  Model  (MEM)  computes  the  population 
distributions  up  to  v  =  40,  however,  the  measurements  extend  no  higher  than  v  =  5  and  more 
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often  only  to  v  =  4.  In  general,  above  some  value  of  v,  the  deactivation  coefficient  yv  is  unity, 
which  occurs  from  as  low  as  v  =  4  to,  perhaps,  v  =  10.  Thus,  it  was  necessary  to  devise  a 
reasonable  procedure  for  extrapolating  the  k,  to  higher  v’s  in  order  to  use  the  MEM  program. 
Two  methods  were  used:  1)  since  the  observed  relaxation  rates  k,  typically  increased  linearly 
with  v,  a  linear  extrapolation  of  k  was  carried  out  until  yv  reached  unity;  2)  the  square  root 
of  the  measured  CARS  signals  were  extrapolated  following  the  nearly-linear  function 
established  by  the  recorded  data.  The  latter  approach  creates  a  renormalization  of  the  data 
such  that  the  recorded  values  ofk  are  shifted  slightly  from  the  values  given  in  Tables  1,  3  or 
5.  In  either  approach,  the  actual  wall  relaxation  rates  k^  are  fitted  to  the  data  from  which 
the  “true”  deactivation  coefficients  yv  are  calculated  using  Eq.(20).  Table  7  gives  the  results 
of  six  studies  of  fitting  the  MEM  to  the  data.  It  should  be  noted  that  these  calculations  were 
very  tedious  and  the  fittings  were  often  plagued  by  ambiguities.  Because  of  the  shear 
complexity  of  the  relationships  between  the  relaxation  rates  of  the  observed  vibrational  states 
and  the  impact  of  the  interplay  between  the  V-V  and  wall  relaxation  rates  on  the  radial 
distribution,  it  was  often  impossible  to  find  a  unique  and  plausible  fit  to  a  given  data  set. 

The  results  in  Table  7  show  a  very  significant  characteristic.  A  surface  giving  high 
wall  relaxation  rates  is  accompanied  by  observed  relaxation  rates  that  are  lower  than  the  wall 

values  This  characteristic  also  carries  over  to  the  deactivation  coefficients  wherein  the  yv  are 

larger  than  the  apparent  values  yva.  This  can  be  seen  for  the  +  column  Pyrex,  the  stainless 
steel  and  the  aluminum  results  given  in  Table  7.  On  the  other  hand,  the  materials  giving  low 
or  weak  wall  effects  show  wall  relaxation  rates  k^,  that  are  very  similar  to  or  smaller  than 
the  observed  relaxation  rates  k  Likewise,  this  characteristic  carries  over  to  the  deactivation 
coefficients.  This  is  the  case  for  neat  Pyrex,  gold  foil  and  Teflon.  These  characteristics 
suggest  that  the  V-V  interactions  have  a  significant  impact  on  the  radial  (i.e.,  transverse) 
distributions  of  the  populations  of  vibrationally-excited  molecules  when  the  wall  interaction 
is  strong.  In  effect,  the  V-V  processes  redistribute  the  populations  so  as  to  make  the 
observed  axial  decay  rates  slower  than  the  wall  rates.  For  weakly-interacting  surfaces,  the 
wall  effects  do  not  create  significant  radial  population  gradients  such  that  the  observed  axial 
decay  of  the  populations  track  closely  the  longitudinal  decay  of  the  populations  at  the  wall. 
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Table  7.  Results  of  Master  Equation  Model  fitting 


Material 

T 

(K) 

P 

CD 

k  or  y 

v=  1 

IH 

IH 

■ 

■ 

Signal 

Extrp. 

Neat 

360 

TT 

Ms'1) 

16±14% 

38±10% 

54±8% 

67±19% 

yes 

Pyrex 

^v(waH) 

24 

43 

53 

77 

to 

Yva 

5.2 

13 

20 

27 

v=9 

Yv 

7.9 

15 

20 

34 

49 

+ 

340 

14 

Ms-1) 

19±1 1% 

38±5% 

57±5% 

75±11% 

no 

column 

^v(waD) 

30 

50 

70 

90 

Pyrex 

104yv, 

6.4 

14.8 

25 

39 

104yv 

11 

21 

35 

55 

99 

Gold 

351 

17 

K  (s'1) 

16.7±5% 

37±3% 

53±3% 

64±3% 

65±5% 

no 

Foil 

^vtwaD) 

25 

36 

48 

58 

104yva 

5.5 

14.5 

24.1 

34 

34 

104yv 

8.9 

13.8 

20.9 

28 

43 

SS  312 

382 

14 

K  (s'1) 

21  ±9% 

50±7% 

73±5% 

111±12% 

yes 

47 

79 

100 

116 

to 

10^ 

8.7 

28 

63 

104 

v=9 

l04yv 

25 

81 

424 

104 

A1  3003 

375 

14 

MM 

8  4±9% 

21  ±8% 

34±7% 

45±8% 

gg!| 

^v(waD) 

17 

33 

46 

63 

WES* 

104yva 

2.9 

7.6 

14 

20 

104yv 

6.1 

13 

21 

34 

54 

■ 

Teflon 

334 

14 

K  (sl) 

5.5±17% 

1 1±19% 

18±14% 

28±17% 

1 

Mvall) 

5.7 

94 

14 

20 

io4yva 

1.7 

3.4 

6.3 

10 

I04yv 

1.8 

3.0 

4.5 

6.8 

8.8 

■ 

XPS  Measurements  on  Selected  Specimens 
Examinations  of  the  surfaces  of  selected  specimens  were  carried  out  by  the  Surface 
Sciences  Laboratory  of  the  University  of  Dayton  Research  Institute.47  The  specimens 
consisted  of  witness  coupons  that  were  inserted  for  a  deactivation  experiment  in  the  catch 
tube  to  intercept  the  afterglow  flow  approximately  3  cm  after  the  laser  probe  volume  or  the 
Pyrex  specimen  tubes  were  broken  to  provide  witness  coupons  for  XPS  (x-ray  photoelectron 
spectroscopy)  examination.  Table  8  lists  the  13  specimens  that  were  examined.  The 
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objective  of  these  measurements  was  to  determine  the  degree  of  modification  of  the  surface 
constituency  by  exposure  to  active  nitrogen  either  in  the  afterglow  or  in  the  positive  column. 
This  concern  was  prompted  by  the  observation  that  the  observed  relaxation  rates  varied  over 
multiple  experiments  on  the  same  Pyrex  specimen  tube. 


Table  8.  Description  of  specimens  measured  with  XPS 


Specimen  No 

Description  of  Specimen 

1 . tube  i 

Pyrex;  inside  near  flow  outlet  end;  Part  A 

2;  tube  i 

Pyrex;  outside  near  flow  inlet  end;  Part  B 

3;  coupon  ii 

Pyrex;  witness  in  afterglow  at  outlet  end  of  test  specimen;  Part  A 

4,  tube  ii 

Pyrex;  positive  column  exposure;  white  deposit  on  outside;  Part  B 

5;  tube  ii 

Pyrex;  positive  column  exposure;  inside  near  flow  outlet  end 

6;  tube  iii 

Pyrex;  long  term  exposure;  inside  near  flow  outlet  end 

7,  coupon  iii 

Pyrex;  to  check  for  metal  deposition  from  discharge  system 

8;  coupon  iv 

Pyrex;  used  with  Au  foil  specimen 

9, coupon  v 

Pyrex,  before  nitriding 

10,  coupon  vi 

Pyrex,  during  nitriding  in  afterglow;  gas  contact  (convex)  side 

1 1 ,  coupon  vi 

Pyrex;  during  nitriding  in  afterglow;  back  (concave)  side 

12,  coupon  vii 

Pyrex;  after  nitriding  in  afterglow 

13;  coupon  viii 

Au  foil,  clipped  from  specimen  at  gas  outlet  end;  side  in  contact 
with  active  nitrogen 

Table  9  summarizes  the  results  of  the  measurements  determined  by  XPS.  This 
technique  detects  the  binding  energy  of  atoms  in  the  surface  region  to  a  depth  of  about  30 
Angstroms  The  recorded  signal  count  rate  for  a  given  energy  is  proportional  to  the  number 
of  bonds.  The  carbon  readings  originated  from  the  material  used  to  mount  the  specimens  and, 
therefore,  were  intrinsic  to  the  XPS  system.  There  was  no  evidence  of  carbon  on  the 
specimens.  A  good  example  of  an  interpretation  of  these  measurements  is  the  study  of  the 
inside  surface  of  specimens  Nos.  1  and  2  in  Table  8  taken  from  Pyrex  specimen  tube  i.  The 
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outlet  in  of  this  tube  had  a  1  atomic  %  of  nitrogen  bonded  to  oxygen  and  1.5  atomic  % 
bonded  to  silicon  while  the  concentration  of  the  inlet  end  was  less  than  0.3  atomic  %.  This 
result  confirmed  our  suspicions  that  the  surface  of  the  specimen  tube  was  being  modified  by 
the  afterglow.  Working  from  the  formula  and  density  of  Pyrex,  the  surface  at  the  outlet  end 
had  a  8.4%  coverage  of  nitrogen  adatoms.  A  number  of  deactivation  experiments  had  been 
run  on  this  specimen  tube  with  an  accumulated  exposure  to  the  afterglow  of  more  than  70 
hours  These  experiments  gave  not  only  different  values  of  observed  relaxation  rates,  k„  and 
the  corresponding  yw  but  different  dependencies  of  yva  on  v.  The  results  from  the  surface 
lab  supports  the  contention  that  the  deposition  of  nitrogen  atoms  on  the  surface  affects  the 
vibrational-  energy-to-surface  exchange  mechanisms.  The  correlation  of  the  XPS  studies  with 
exposure  to  active  nitrogen  points  to  two  observations:  1)  continued  exposure  tends  to 
increase  yla  and  2)  the  increase  of  yva  with  v  relative  to  yla  decreases.  As  a  general  rule,  the 
observed  deactivation  rates  k,  increased  linearly  with  v;  however,  as  the  time  of  exposure  to 
active  nitrogen  increased,  the  slope  of  k,  with  v  decreased. 


Table  9.  Atomic  percent  of  surface  composition  of  specimens  given  in  Table  8 


No. 

c 

0 

Si 

N-O 

N-Si 

Na 

Cu 

Ni 

Fe 

A1 

F 

Au 

1 

12.1 

Kia 

warn 

warn 

2 

14.3 

60.1 

25.1 

<0.3 

<0.3 

3 

19.2 

51.0 

19.5 

— 

,0.5 

— 

5.8 

5.0 

4 

36  6 

38.5 

8.2 

— 

0.7 

15.7 

<0.3 

5 

7.4 

62.9 

27.3 

1.1 

1.3 

6 

13.8 

60.5 

24.3 

1.2 

<0.3 

7 

21.9 

55.1 

18.9 

— 

<0.5 

1.0 

2.5 

8 

26.1 

46.2 

17.4 

— 

— 

5.8 

2.3 

— 

— 

— 

— 

2.2 

9 

14.7 

58.6 

22.1 

— 

<0.6 

— 

— 

1.1 

<0.2 

msm 

— 

0.3 

10 

11.5 

52.3 

26.1 

1.8 

ma 

— 

— 

<0.2 

<0.3 

<0.5 

— 

— 

11 

27.1 

44.9 

27.0 

!  — 

<0.6 

14.8 

— 

<0.2 

— 

0.8 

— 

— 

12 

53.5 

2.0 

1.0 

HI 

— 

— 

— 

0.8 

irai 

Et! 1 

lUrl 

... 

Belli 

— 

EiKfl 

■ran 
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CONCLUSIONS 


The  results  presented  in  this  work  represent  the  first  extensive  study  of  deactivation 
of  vibrationally-excited  molecules  for  the  v  =  1,  2,  3,  4  and  5  hot  bands  on  a  number  of 
different  surfaces.  Moreover,  this  work  includes  the  first  attempt  to  fit  a  relatively  complete 
set  of  rate  equations  for  the  prevalent  processes  to  such  data  in  order  to  extract  the  true 
surface  deactivation  coefficients.  The  values  reported  here  for  Pyrex  are  consistently  lower 
for  the  v  =  1  band  by  a  factor  of  about  two  than  the  values  published  by  Black  et  al.8  which 
can,  at  least  in  part,  be  attributed  to  their  lower  temperature  of  300  K.  This  is  supported  by 
the  temperature  behavior  of  the  results  obtained  with  the  gold  foil  wherein  the  yla  decreases 
by  almost  a  factor  of  four  for  a  change  in  temperature  of  56  K.  The  results  obtained  for 
metals  were  qualitatively  the  same  as  reported  by  Black  et  al.  although  our  values  for 
aluminum  were  somewhat  lower,  which  could  be  due  to  different  alloy  used  here  and  our 
higher  temperatures  There  is  a  considerable  discrepancy  between  their  value  for  Teflon  and 
ours  With  the  exception  of  titanium,  Teflon  gave  the  lowest  yla  values  while  Black  et  al. 
found  a  value  for  Teflon  equal  to  that  for  Pyrex.  The  results  obtained  here  for  titanium  were 
quite  surprising  and  represent  a  significant  new  result.  The  obvious  implication  is  that  a 
container  of  vibrationally-excited  nitrogen  made  out  of  titanium  could  provide  some 
significant  degree  of  storage  of  this  energy. 
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